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V.l HEIGHTS 


There are several subroutines in GASP that srs used to determine 
the ccnponsnt weights of the aircraft configuration. Those can bo grouped 
into the propulsion related conponants end the airframe components. The 
airframe group also includes subroutines that deal with design load conditions* 
aircraft balance* and tail sising. Figure V.1.1 liets the weight related 
subroutines. 

There are several options available for determining propulsion related 
weights and the options depend on whether it is a turbine or internal combustion 
type of engine. All turbine engines' weights are determined in subroutine 
ENGWGT. internal combustion type engine weights may be determined from 
trend equations in ENGWGT or from a more detailed formulation as found in 
HOPWSZ and RCWSZ subroutines . Propeller weight may be determined by ENGWGT 
or the weight methodology of Hamilton-Standard which is contained in subrou- 
tine WAIT. If a known propulsion system is being studied* the component 
weights of the propulsion system may be input directly. 

The weight of each aircraft structural component is estimated by a 
regression analysis equation derived from similar aircraft* and a large 
number of parameters is used to describe the aircraft. These include linear 
dimensions, dimensionless parameters such as aspect ratio and fineness ratio* 
and cargo data such as number of passengers. Their use is based on the 
assumption that aircraft component weights tend to be in proportion to certain 
of its dimensions. The tail sizing subroutine is an option which develops 
horizontal and vertical tail dimensions which satisfy specified stability and 
control criteria. 
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PROPULSION; 

ENGWGT (160) 

CONTROLS AND COMPUTES ENGINE WEIGHT 
COMPUTATIONS 

HOPWSZ 

HORIZONTAL OPPOSED PISTON ENGINE WEIGHT 


AND SIZE 

RCWSZ 

ROTARY COMBUSTION ENGINE WEIGHT AND SIZE 

WAIT 

PROPELLER WEIGHT 

AlfiEBAHE: 

WGHT (510) 

STRUCTURAL WEIGHT AND WING POSITIONING 

DLOAD 

DESIGN LOAD FACTORS 

TAIL 

TAIL SIZING FOR LONGITUDINAL AND DIREC" 


TIONAL STABILITY 


FIGURE V.1.1 WEIGHT SUBROUTINES 
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V.1.1 Propulsion System Heights 

V.l.1.1 Subroutine ENGWGT 

This 160 card subroutine is eoncemsd with the weight of the propulsion 
system, and scow 30 input parameters are used to calculate any of about 20 
output quantities. The steps followed depend strongly upon the type of engine 
and propeller being used, and frequently the weights are expressed as statis- 
tical functions of other parameters of the aircraft. 

The call integers of the subroutine are NTYEX and NTYPX which correspond 
to engine and propeller types as follows i 


NTYEX - 


1, 2, 3, 11, 12, 13i reciprocating engines 

4, 14 t rotating combustion engines 

5, 6 i turboshaft and turboprop engines 

7: turbojet or turbofan engine (V.1.1) 


NTYPX - 



lit 

12t 

13 1 
14t 


5, ISt* 

6, lot 


fixed-pitch propeller 

constant speed propeller 

constant speed, full-feathering propeller 

constant speed, full- feathering, deice 

propeller 

constant speed, full feathering, reverse 
deice propeller 

Q-fan propulsor (V.1.2) 


The subroutine's generality is possible only because of the large number of 
inner paths depending on these two integers, and it is not feasible to describe 
in detail the flow of logic for every pair NTYEX and NTYPX. 
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The intagtri NTYEX ■ 11 to 14 provide tho option to un subroutines 
HOPffSZ and FCWSZ to dovolop par am* trie valuta of weight and also of horizontally 
oppoaad piaton anginaa or rotary combustion type of anginaa . Similarly, valuta 
of NTYPX greater than 10 utilize aubroutine WAIT through a call to ENGINE to 
return valuer of propeller weighta. Otherwise, the flow of logic ia parallel 
for arbitrary valuta of the call integers. 

Following the initial zero epecification of component weighta, the 
aubroutine defines a value for sea level static engine specific weight. If 
this has not been input, it is computed by 


Stf 


SLS 


1.5(1 ♦ .15 K spcHG > 

1 * ,2 k spchg 
.5 

.13 


NTYEX » 1, 2, 3 
NTYEX - 4 
NTYEX - 5, 6 

NTYEX - 7 (V.1.3) 


This parameter has units of lb per HP except for the turbojet or turbofan 
engines, NTYEX ■ 7, in which case the units are lb per lb of thrust. The 
same parameter is returned by subroutines HOPWSZ (NTYEX ■ 11, 12, 13) and 
RCWSZ (NTYEX - 14), or is equal to .5 (NTYEX - 15, 16). 

In terms of this parameter, the engine weight is 


W ENG " SW SLS ^MSLS 


(V.1.4) 


or 


W ENG “ SW SLS ^SLS 


(V.1.5) 


where 

HP MSLS " maximum sea level static horsepower 
FN„. _ ■ sea level static thrust 
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Tha nacslls Might is than developed ss ths product of two paramo torn 


V NAC " ^NAC '\*AC 

whsro 

uw NAC - nacelle weight per unit surface area, 
lb per sq ft 

Vc * tot * 1 n-.elle area, sq ft 


Zn this equation UW fjAC is an input parameter and A^, is computed internally 
(KNAC - 0, 1) or input (KNAC - 2). 

Next, if the number of engines in two or more (EN p f 1) , the pylon weight 
is expressed as 


where 


.736 

“ F___ ♦ W ) 

PYLON PYL ENG NAC 

F_„, ■ input scale factor 
PYL 


(V.1.7) 


and W ' W are obtained from equations (V.1.4) or (V.1.5) and (V.1.6). 

ENG NAC 

For propeller types NTYP - 1 to 6, the propeller weight is given by 
numerically complex functions of propeller diameter °p R Qp • numJaer of 
blades b l# propeller blade activity factor A^,, engine speed RPM and horse- 
power HP , and cruise Mach number EM__„ . These input parameters are used 
to define numerical values 'for the constants XK^, XC W » F^, F^,...., F & . For 
NTYP • 1 to 5, the propeller weight is then estimated by the function 


and for NTYP 


„ 2 .7 .75 15 .12 .5 A 

“fan • SK 1 F 1 F 2 P 3 P 4 P 5 P 6 * XC W 

6, the fan weight is similarly written as 


w, 


FAN 


1.85 .7 .6 .5 .12 .5 

X K1 F 1 F 2 P 3 F 4 P 5 * 6 


(V.1.8) 


(v.i.y) 
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and tha shroud weight is estimated by 

m sh» - *“2 “prop ,V - 1 -* 0 ' 

wh “* “2 u * ““ rlc * 1 f ““ r *" d “prop “ ‘ ,rop * 11 “ dl *“ t “- *»• 

Torque is computed by its definition in terms of horsepower and RPM, 

TORQUE - 550 HP^^/C 2 - II*RPM/60.> (V.1.11) 

Prom this value the gear box weight is estimated by 

WT n - XK, TORQUE* 84 (V.1.12) 

Go J 

where XK 3 ■ .085# unless input as nonzero. Zf the Hamilton-Standard option is 
specified, subroutine ENGINE returns values for the weight parameters which 
are computed in subroutine WAIT. 

For turbojet or turbofan engines. NTYE ■ 7, and the engine weight 
parameters are very simply given between statements 500 and 600 in terms of 
input parameters . i . e . , 


FN„, _ ■ WA SFN_ r _ 
SLS SLS SLS 


* SW„ e FN Te , 
ENG SLS SLS 


nacelle weight. 


Vc ‘ °Vc < n d e “bde* • 


pylon weight. 


W_ 


p /u + u ) 
PYL ' ENG NAC 7 


.736 


PYLON 10 (if single engine configuration) 


(V.1.13) 

(V.1.14) 


(V.1.15) 


(V.1.16) 
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Nhtfn Input piraMteri art 

WA_. „ ■ sea level static airflow 

SLo 

SPN - sea level static specific thrust 
SLS 

SW_ T _ - sea level static engine specific weight 
SLS 

FN__ - sea level static thrust 
SLS 

UW„. _ - nacelle weight per unit area 

NAC 

P • factor for pylon weight 

The remaining statements of the subroutine are concerned only with print 
statements for several values of NTYE and NTYP. 

V.l.1.2 Subroutine HOPWSZ 

This subroutine deals with numerical weight and size of horizontally 
opposed piston engines. The input to the subroutine is descriptive information 
as to engine geometry of which the following are the most important: 

SKWGT - engine weight calibration factor (nominally 1) 

SKWDTH • engine width calibration factor (nominally 1) 

KWH * ratio of engine width to height (nominally 1.3) 

HPMSLS - maximum sea level standard horsepower of engine 
NCYL - number of engine cylinders where 4 $ NCYL $ 8. 

HPQAB ■ ratio of rated horsepower to piston bore area, hp per sq in 
Several numerical scale factors are defined as functions of NCYL and other 
input parameters, after which the principal output quantities are found as 
SWSLS - specific weight (sea level standard) lb per hp 
ANAC - nacelle surface area, sq ft 
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XLN ■ length of nacelle, ft ■ 1.2 * length of engine (XL) 

WN ■ width of nacelle ■ 1.08 * width of engine (W) 

HN - height of nacelle ■ 1.2 * height of engine (H) 

Using these definitions the engine specific weight is given by 


» 


SWSLS - SKWGT * SWN 


( 1800 ^ ) ^ 
1 PSTSPD BMEP 7 


(V.1.17) 


wheze SWN is found by a table look-up. The variable PSTSPD , piston speed, 
is computed from 

PSTSPD - 33000. * HPQAB/BMEP (V.1.18) 


Engine width, W, is given by 


W - SKWDTH * (k * HPMSLS + 30.0) 


(V.1.19) 


k 


.0167? NCYL - 4 
.0125 NCYL - 6 
.0106 NCYL - 8 


(V.1.20) 


For unsupercharged engines the engine width to length ratio, KWL, is computed 
from 

RWL - B - .09 * RWH (V.1.21) 

where depending on engine type 


1.31 or 1.11; 

NCYL - 4 


1.10 or 0.92? 

NCYL - 6 


1.0 or 0.83; 

NCYL - 8 

(V.1.22) 
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For supercharged engines 


RWL - B - .156 * RWH (V.1.23) 

where 



1.53i 

NCYL ■ 4 


B - 

1.40; 

NCYL - 6 



1.31; 

NCVL - 8 

(V. 1. 2d ) 

Engine length and height are 

then given by 



L 

- W/RWL 

(V.1.25) 

and 





H 

- W/RWH 

(V.1.26) 

Nacelle surface area is 

computed from 


AN AC • 

- 0.8 

* <WN + 2.76 * HN) * XLN/144. 

(V.1.27) 


V.l.1.3 Subroutine RCWSZ 

This subroutine deals with weight and sire characteristics of the rotating 
combustion engine. The subroutine begins with a tabulation of several input 
data parameters, after which the input/output computations begin, in much tin- 
same order as was described for the previous subroutine. Principal input 
parameters are listed in the call statement as 

HPMSLS - maximum sea level static horsepower 

ROTN - number of rotors 

GR - ratio of propeller rpm to maximum engine rpm 

SKWGT - engine weight calibration factor (nominally 1.) 
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SKDXAM ■ engine diameter calibration factor (nosdnally 1.) 

SWSLS - specific weight, standard sea level, lb per hp 

Other input quantities passed through COMMON statements, deal with the 
technology date level (IDATE - 1970 or 1980) , and the absence or presence 
of supercharger (KSPCHG ■ 0 or 1) . 

The computational results are straightforward and vary principally 
with horsepower and number of rotors. Output quantities of the subroutine 
are 


GR - ratio of prop rpm t,; engine rpm (operating condition) 

SWSLS ■* specific weight with supercharger, lb per hp 
ANAC ■ nacelle wetted area, sq ft 
XLN - nacelle length, ft. 

V. 1.1.4 Propeller Weights, Subroutine WAIT 

Propeller weight is estimated in this 30->card subroutine as a numerical 
function of seven input parameters: 

NTYP - IWTCON - airplane propeller type (1 to 5) 

ZMWT ■ Design cruise Mach number 
BHP - brake horsepower 
DIA - propeller diameter, ft 
AFT ■ activity factor per blade 
BLADT - number of blades 
TIPSPD- tip speed, ft per sec. 

Then, according to several straightforward but nonlinear functions, the output 


parameters are simply: 
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WT70 • propeller weight, 1970 technology, lb 
WT80 - propeller weight, 1980 technology, lb 


Equations employed are 



where 


(V.1.28) 


WT70 or WT80 ■ W T ■ propeller wet weight, lbs. (excludes spinner, de- 
icing and governor) 

DZA • D » propeller diameter, ft 
FLADT • B * number of blades 
AFT ■ A.F. ■ blade activity factor 

Vmjp 

N ■ propeller speed, rpm (take-off ■ ■ - ■ ; V ■ TZPSPO) 

WQ TIP 

BHP » SHP ■ shaft horsepower, HP (take-off) 

ZMWT - M - Mach number (design conditions maximum power cruise) 


V 


C W* y (lo) ( 8 )(im) ( i2 »5 22 ) ■ Counterweight wt., lbs. 

(V.1.2^) 

u, v, and y values for use in the weight equation are taken from 


the table below 



Propeller 


1 

Type 

Technology | 

(NTYP) 

1970 

1980 

1 

(1) 

(1) 

2 

(2) 

(2) 

3 

(3) 

(3) 

4 

(3) 

(4) 

5 

(3) 

(5) 



** 

u 

V 

19 

(1) 

170 

m 

0.35 

0 

(2) 

200 

Q 

0.35 

0 

(3) 

220 


0.40 

5.0 

(4) 

190 


0.40 

3.5 

(5) 

190 

Hi 

0.30 

0 


Propeller types associated with above , u, v, y are as follows: 

1. all fixed-pitch props 

2. McCauley non-counterveighted , non- feathering, constant speed props 

3. all Hartzell, all Hamilton-Standard small props, and feathering 
McCauley 

4. fiberglass-bladed, constant speed, counterweighted, full feathered 

5. fiberglass-bladed, constant s.oeed, double-acting (non-counterweighted) , 
full feathered, reverse 

V.1.2 Airframe Weight 

The weight of the structure, flight controls, and payload components 
of the aircraft are determined in subroutines WGHT (Section V.1.2. 2). 

Subroutine DLOAD (Section V. 1.2.1) determines the minimum design speeds, and 
structural load factors while subroutine TAIL (Section V.1.2. 4) is used in 
one of the stability and control options. 
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V. 1.2.1 Subrouting DLOAD - Design Speeds and Load Factors . The 


purpose ox this subroutine is the computation of minimum design airspeeds 
and load factors of an aircraft in accordance with Federal Aviation Regulations , 
(FAR) Part 23 or 25. These speeds and load factors jure used to permit 
computation of structural weights by subroutine WGHT. Both maneuver load 
factor and gust load factor are considered. 

The primary input parameter to this subroutine is the airplane struc- 
tural design category indicator CAT^, defined as 


cat d 


0, FAR Part 23 normal category 

1, FAR Part 23 utility category 

2, FAR Part 23 aerobatic category 

3, FAR Part 25 transport category 


If FAR Part 23 requirements jure selected, wing loading (WOS) is also an 
importjutt parjuneter for determining the minimum design airspeeds. Other input 
parameters passed through the subroutine argument list lure 


V MLFSL " Est ^ inate structural design equivalent airspeed, mph 

C “ curve slope, per rad. 

LALPH 

C_ “ wing mean aerodynamic chord, ft. 

BAKW , 

The subroutine begins by changing the units of the estimated structural 
design velocity from miles per hour to knots, 

V -V /1.15 

MLFKT MLFSI/ 
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and if CAT d - 3 branches to statement 50 for transport category aircraft. 
If one of the PAR Part 23 categories has been selected* the minimum design 
airspeeds are determined and checked. 


For CAT d ■ 0 or 1 i the minimum design cruise speed is 


CMIN 


33. ^toT, WOS < 20 

[33 - .055 (WOS- 20) ] ^OS WOS > 20 (V.1.30) 


and for CAT Q ■ 2, 


CMIN 


36 VWOS WOS < 20 

[36 - .0925 (WOS- 20)) ^TOS WOS > 20 (V.1.31) 


The minimum design dive speeds are also functions of the wing loading; 
i . 6. * 

Normal Category, CAT Q ■ 0; 


DMIN 


1.4 V 


CMIN ' 

[1.4 - .000625 (W0S-20))V 


CMIN 


WOS < 20 

WOS 2 20 (V.1.32) 


Utility category, CAT Q ■ Is 

1.5 V 


CMIN 

'DMIN f [It 5 - .001875 (WOS-20))V 


WOS < 20 


CMIN WOS > 20 (V.1.33) 


Aerobatic category, CAT Q - 2: 


DMIN 


1.55 V 


CMIN' 


WOS < 20 


[1.55 - .0025 (WOS- 20) JV _ WOS > 20 (V.1.34) 

CHIN * 


V-l 

14 


For transport category (CAT^ ■ 3) thara are no lower limits established 
for the design cruise speed as for FAR Part 23 aircraft. The minimum design 
dive speed for transport category aircraft is 

V DMXN " 1,2 V MLFKT " 1,2 V M0 (V.1.35) 

The limit maneuvering load factors pex the FAR's are then specified as 

3.0, CAT d - 0 

4.4, CAT d - 1 

6.0, CAT d - 2 

2.5, CAT q - 3 (V.1.36) 

The gust load is evaluated at the altitude at which maximum operating equivalent 
airspeed (VMO) is equal to the speed for maximum operating Mach number (M Mq ) 
so long as the altitude falls in the band 



0 ✓ h < 20,000 feet (SIGMA > .53281) (V.1.37) 

v V 


For a non- transport category aircraft (CATD ft 3) this limit is 12,500 feet. 
(SIGMA > .6020). 

The gust load factor is calculated for a 50 ft/sec gust at VCMIN and 
a 25 ft/sec gust at VDMXN, and the most critical condition selected for the 
gust load condition. 


max (G, 


g lfo } d 


(V.1.38) 


where the design cruise speed and design dive speed load factors are 


G LFC ■ 1 * 50 g lfk c lalph v cmim /< 490 - mos > (v - 1 - 3,) 

G LFD 1 + 25 G LFK C LALPH V DMIN /(498 * W0S ^ (V.1.40) 
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Xn t he** equation* , 

- .88 Ug/tS.X ♦ U G ) (V.1.41) 

where the airplane mass ratio is defined as 

U Q - 2. WOS/(.00237*SIGMA*C bark ^^^^*32. 2) (V.1.42) 

The closing computation deals with the ultimate load factor, which is 
simply 50 per cent greater than the moat critical of the maneuver or gust 
load factor 

_ 1.5 EM,, . 

" 15 V • Otp a 1*1! 

V.l.2.2 Subroutine WGHT . This subroutine is used to compute the weight 
of the structural and flight controls components of the aircraft using the 
geometry established by subroutine SIZE. It is also possible to size tip 
tanks if the wing volume is inadequate for carrying the available fuel. An 
additional option provides for relocating the wing on the fuselage and 
resizing the horizontal tail to provide a given level of longitudinal stability. 

The weights of aircraft structural and flight controls components arc 
estimated using nonlinear statistical weight equations, and the inner loop 
is used to compute fuel weight, which may require the consideration of wing 
tip tanks, as shown in Figure V.1.2. A secondary loop involves the resizing 
of the engines, and the third major loop requires that the aircraft center 
of gravity location satisfy certain stability conditions. This may mean that 
the tail moment arm must be changed which requires that the inner loop be 
recomputed. 

V-l 
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Input 

Gross Weight 
Payload 
Geometry 
Design Speeds 


ompute Weights 



r Resize 
Engines 


/ Tip ^ 

Tank Fuel 
s£onvcrge« 


Engines 

Resized? 


Recompute 

Drag 



Balance Aircraft: 
Recompute Tail Are 



— Tai 1 \ 
Areas Con- 
vyerge? 



Output 

Component Weights 

Fuel 

Loads 

Wing Locat 



Figure V.J,2 Simplified Flo*/ Chart for Subroutine WGHT 
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The complexity of the problem is increased by tha fact that wing weight 
depanda on fuel weight, and becauae changaa may require the conaideration of 
wing tip tanka, aa shown in Figure V.1.2. A secondary loop involves the 
resizing of the engines, and the third major loop (IF LOWING f 0) requires 
that the aircraft center of gravity location satisfy certain stability 
conditions. This may mean that the tail monenc arm must be changed which 
requires that the inner loop be recomputed. 

The complexity of vhe problem is further increased by the fact that wing 
weight depends on fuel weight, and because changes in the wing tip tanks 
may require the engines to be resized. This nonlinear dependence of the 
aircraft parameters also shows in the balance option, for which a change in 
tail area and moment arm also varies the weights of both wing, tails and fuel. 
As suggested by Figure V.1.2 the following four headings will be discussed 
here: 

1. Component weights, 

2. Fuel weight, 

3. Tip tank weight, and 

4. Center of gravity 

The subroutine begins by defining a number of numerical parameters, 
including fuel density FUEL^ and quarter chord sweep angles of wing and tail 
surfaces. The lift curve slope of the wing is 

C LALPH ■ n * AR/(1. + 1 *vC^r 2 > (V.1.44) 

where C^ gnd C 2 are functions of quarter chord sweep, aspect ratio, and 
cruise Mach number: 


V-l 

17A 



(V.1.4S) 


C x - lAR/(2. CO»(R U|c4 ))) 2 

c j - 1- - 1“ C RU CO,(R Ue« ,lJ (V.1.46) 

The subroutine DLOAD is then called, which returns values for the last 
six call arguments 

V dmxn " design dive speed, kts 

Ojjp ■ ultimate load factor, g's 

- gust load factor, g's 
EM^p ■ maneuver load factor, g's 

- maximum allowable operating airspeed, kts 
EM^ ■ maximum allowable operating Mach number 

The input quantities to DLOAD are the first five arguments, which are 
maximum velocity VM^^, wing loading W , FAR structural design category 
CATp, lift curve slope C LALpH , and wing mac, C BARW ' itaa COMMON/SIZE/ . 


V. 1.2. 2.1 Fixed Equipment Weight. Unless input in namelist INGASP 
as WFEX, the weight of fixed equipment is estimated in pounds as 

W„„ ■ 61.75 * SEATS 2 - 352.5 * SEATS + 533. (V.1.47) 

FE 

where 


SEATS ■ PAX + 1 (V. 1.48) 

and PAX is the number of passengers carried by the- aircraft. This trend is 

applicable for two-seat trainer configurations through 8 to 10 place business 

type aircraft. Figure V.1.3 shows the items that are considered a part of 

the fixed equipment as well as the breakdown for some current aircraft. 

Program modified Feb. 1980 replacing equation 
V-l V.1.A7. If WFEX not input WFE and WFUL (next 
10 section, pg. V-l-20) are computed by new sub- 
routine WFIXEU. Figures V.1.3 and V.1.4 are 
still useful guidelines. 
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AIRCRAFT 




ITEM 

Cessna 

150 

Piper 

Arrow 

Cessna 

210 

Cessna 

340 

Citation 


DHC-6 


Auxilary Propulsion Unit 

- 

- 

- 

- 

- 

- 

- 

345.5 

Instrument & Navigation 

5.8 

15.4 

20.0 

61.9 

85.7 

80.0 

70.4 

211.7 

Hydraulic and Pneumatic 

3.3 

19.8 

51.1 

15.3 

94.2 

106.0 

43.3 

236.7 

Electrical 

44.7 

49.6 

59.5 

131.5 

321.4 

441.0 

343.4 

1073.7 

Avionics 

- 

- 

- 

1.1 

361.5 

300.0 

65.0 

607.1 

Furnishings 

45.1 

117.4 

195.0 

244.9 

688.0 

488.0 

785.7 

1981.9 

Air Conditioning 

5.4 

- 

7.0 

142.1 

187.8 

200.0 

243.3 

497.2 

Anti- Ice 

- 

- 

2.9 

3.7 

63.0 


190.2 

284.2 

Auxilary Items 

- 

- 

2.5 

2.7 

2.4 

- 

- 

6.0 

Paint 

- 

13.3 

13.6 

37.5 

36.0 

50.0 

38.2 

5.9 

TOTAL (W ) 
FE 

104.3 

215.5 

351.6 

640.7 

1840.0 

1665.0 

1779.5 

5249.9 


FIGURE V.l 3 WEIGHT OF FIXED EQUIPMENT 












V. 1.2. 2. 2 Fixed Useful Load. The weight of fixed useful load. WFUL, 
must be input in namelist INGASP. Figure V.1.4 shows the items that are 
usually considered a part of this category as well as the breakdown for some 
current aircraft. 


V. 1.2. 2. 3 Payload Weight. The program is set up to determine the 

range capability of the configuration for three payloads. The three conditions 

are known as the maximum payload, maximum fuel, and the design payload. These 

are outlined in the following paragraphs. 

The maximum payload is determined from the number of passengers (PAX) 

and unit weight of passenger and luggage (UW ) as input in namelist ingasp; 

PAX 

i . 6 . , 


V* - PAX (V. 1.49) 

PL PAX 

where 0W pw ■ 200 pounds is the typical default value. 

The maximum fuel case payload is determined by the total fuel volume 
capacity of the configuration, constrained by the gross weight or usable- 
fuel volume of the aircraft. The mathematical representation is 


W BI - H 

PL G 0E FMAX 


(V.1.50) 


where w„ - design gross weight, lb 
G 

W QE ■ aircraft o^ -ating weight empty, lbs. 

■ maximum fuel weight, lbs. 

FMAX 

If W pL is negative, the payload would be set to zero and W FMAX reduced to 
stay below the gross weight limit. 

The last case is the design payload case. The design payload weight is 
input as W in namelist INGASP. If no value is input, then the maximum 

PLX 
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AIRCRAFT 

ITEM 

Cessna 

150 

Piper 

Arrow 

Cessna 

210 

Cessna 

340 

Citation 

HffHB 

DHC-6 

Gulfstream 

I 

Crew & Crew Baggage 

200.0 

200.0 

200.0 

200.0 

340.0 

385.0 

400.0 

420.0 

Trapped/Unusable Liquids 

- 

- 

6.0 

41.0 

109.9 

1j2.0 

35.0 

229.0 

Engine 0.1 

*11.3 

15.0 

22.5 

49.0 

37.0 

15.0 

54.0 

113.0 

Passenger Service Items 

- 

- 

- 

- 

- 

10.0 

- 

145.0 

Optional Equipment 

— 

30.9 


247.0 

155.3 


| 


Total (W^) 

211.3 

245.9 

228.5 

537.0 

652.2 

542.0 

489.0 

907.0 


FIGURE V.1.4 HEIGHT OF FIXED USEFUL LOAD 






















payload case previously deacrlbad becomes the design condition. 

V. 1.2. 2. 4 Propulsion System installation. Following statement 3, 
the propulsion system weights are accounted for as 


W ■ primary engine weight 

Er 


W„„_ F.N 
ENG P 


W PROP " P r °P UlB ° r W6ight " EN P W PRO Pl 
w ■ primary engine section structural weight 

• Eb 

<Vc * “pylon 1 ™, 

W PEI “ P rimar y engine installation weight • s Kp EI W gp 

WP„ m - CK_ + CK_ W_ + W 

STAR 5 EP 7 PEI PROP 


(V.1.51) 

(V.1.52) 

(V. 1 . 53) 
(V.1.54) 
(V.1.55) 


where w PROP^ » W NAC * W pyjx)N have been determined by subroutine ENGWGT 

or input. SK p EI is input in namelist INGASP to account for engine installation 

hardware. It may be input as zero if it has been accounted for elsewhere. 

CKg and CK^ are set to 1.0 internally-. 

If the engine section weight W does not have a value at this point, 

PES 

it will be determined by 


W PES " SK PES W EP 


(V.1.56) 


where SK p ES i® input as .338, typically. 

V. 1.2. 2. 5 Landing Gepr. For the normal tricycle gear geometry, the 
total landing gear weight including the running gear (wheels, tires, brakes, 
etc.), structure (shock struts, drag struts, support structure, etc.), 
and controls (retraction, steering, systems, etc.) is expressed as a percentage 
of the design gross weight where 


W LG • (Slt L0 ,W 0 


(V. 1.57) 
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where w - total weight of the landing gear 

SL m landing gear 
^LG gross weight 

W - design gross weight 
G 


(V.1.58) 


The percentage will vary between 0.015 and 0.080 depending on the complexity 

and design loads of the system. Conventional landing gear with retracting 

systems, operating on improved runways, normally run between 0.025 and o.o ' 

STOL-type systems operating on rough runways require longer and larger 

alighting gear components to accomodate the aircraft's higher rotational 

angle and sink speeds required to operate at the shorter field lengths. 

SK for the STOL aircraft will normally vary between 0.035 and 0.08. The 
LG 

main gear usually weighs about 80 per cent of the total gear weight. The 
SK term in the weight expression above is the value that is input in 
namelist INGASP or a default value of 0.0318 is used. The weight fraction 
of the main gear is determined by inputting SK^ G or using the default value 
of 0.80. 

Figures V.1.5 to V.1.7 are included as a guide in selecting SK^. It 
includes the total gear weight as a fraction of the gross weight for a 
sampling of military, commercial, and general aviation aircraft. 

V. 1.2. 2. 6 Empennage Weight. Various geometric parameters for the 
vertical and horizontal tail then follow. The tangent of the leading edge 
sweep angle of the vertical tail is 


TSWVLE - 


1-sus, 

u+suyAR^ 


+ tanOWPQCy) 


where SLM^ 

ar vt 

SWPQCy 


■ taper ratio of vertical tail 

• aspect ratio of vertical tail 

* quarter chord sweep of vertical tail 


(V.1.59) 
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Aircraft 

Gross Weight 
(lb.) 

Total 

Gear Weight 
(lb.) 

Fraction of 
Gross Weight 

iSK tG> 

Bede SJ 

1,000 

41 

0.041 

Cessna 150 

1,600 

114 

0.07125 

Cessna '72 

2,300 

117 

0.05087 

Piper Arrow 

2,600 

98 

0.0377 

Cessna 182 

2,800 

134 

0.0478 

Cessna 210 

3,800 

182 

0.0479 

Cessna 340 

5,975 

268 

0.0448 

Nomad 22* 

8,000 

400 

0.05 

Citation 

11,650 

425 

0.0365 

DHC-6-330* 

12,500 

606 

0.0485 

Merlin IV 

12,500 

623 

0.0498 

Learjet 24 

13,500 

429 

0.0318 

Fregate 

23,810 

1,113 

0.0467 

Gulf stream I 

35,100 

1,237 

0.0352 

Jetstar I 

40,921 

1,081 

0.0264 


FIGURE V.1.5 GENERAL AVIATION LANDING GEAR WEIGHTS 


•STOL-type aircraft 
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Aircraft 

Gross Weight 
(lb.) 

Total 

Gear Weight 
(lb.) 

Fraction of 
Gross Weight 

isk lg' 




F-27 

35,700 

1,884 

0.0528 

VFW614 

44,000 

1,620 

0.0368 

Convair 440 

49,100 

2,158 

0.0439 

F-28 

65,000 

2,649 

0.0407 

737-200 

100,000 

4,038 

0.0404 

DC-9-32 

108,000 

4,182 

0.0387 

727-100 

161,000 

7,211 

0.0448 

Convair 880 

184,500 

6,933 

0.03758 

DC-8-62 

335,000 

11,449 

0.03417 

707-320 

336,000 

12,982 

0.0386 

DC- 10- 10 

430,000 

18,581 

0.0432 

747 

775,000 

32,220 

0.0416 


FIGURE V. 1.6 COMMERCIAL LANDING GEAR WEIGHTS 
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Aircraft 

Gross Weight 
(lb.) 

Total 

Gear Weight 
(lb.) 

Fraction of 
Gross Weight 

<sk lo> 




DHC-4* 

28,500 

1,398 

0.049 

Breguet 941* 

46,000 

2,626 

0.0571 

DHC-5* 

41,000 

1,800 

0.0439 

C-123B 

54,000 

2,331 

0.0432 

C-118 

107,000 

3,895 

0.0364 

C-130B* 

135,000 

4,824 

0.0357 

C-130E* 

155,000 

5,077 

0.0327 

C-124C 

185,000 

11,700 

0.0632 

C-133A 

275,000 

10,635 

0.0387 

C-14IA 

316,000 

10,529 

0.0333 

C-5A 

728,000 

37,628 

0.0517 


FIGURE V.1.7 MILITARY AIRCRAFT LANDING GEAR WEIGHTS 


*STOL-type aircraft 
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The vertical tail root chord and horizontal tail root chord ia then expressed 
as a function of area, span, and taper ratio. 


CRCI^ - 2 S^tB^d + suy) (V. 1.60) 

crcl ht “ 2 s ht /ib ht (1 4 “Si 11 


where 

B ■ span of vertical tail 
VT 

“ the planform area of vertical tail 

Bj|t ■ span of horizontal tail 
S - the planform area of horizontal tail 

HT 

Similar expressions follow for the chordwise distance from root chord 
leading edge to mean aerodynamic chord leading edge, XV^^, and for the 
distance from aircraft nose to wing center of gravity, which is shown as 

el wing in Figure v.i.s 


The weights of the horizontal and vertical tails are determined from 
the weight trend equations developed in Reference 1, and presented below. 
Horizontal Tail 

W HT “ 350 00 054 (V.1.61) 


f h s ht log( vd min ) 

100*ELT h TC ht tRCL HT 


(V.1.62) 


and 


F H - 10 ‘ 6 W G sk y b ht sk tl (1 + 2 * SL V/(1 + slm h ) 


(V.1.63) 


Vertical Tail 


- 380 (k) 


0.54 


where 


k - 


(F y ♦ S 


AH V 2> S VT 


log (vd min > 


100 ELT„ TC CRCL 

V VT VT 


(V.1.64) 

(V.1.65) 
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and 


-6 


r y - .5 X 10 W G SK z (Ely ♦ B) ^ (1 ♦ 2 « SIK y )/(l ♦ Sl/y 

(V.1.66) 


where the following symbols have been used i 

S ■ horizontal tail planform area, sq. ft. 

HT 

S VT - vertical tail planform area, sq. ft. 

» VD U _„ - dive velocity, knots 

MIN 

S - location of horizontal tail on vertical tail (fraction of 

AH 

vertical tail span from root of vertical tail) 

ELT - horizontal tail moment arm, ft. 

H 

ELT V - vertical tail moment arm, ft. 

TC ht ■ horizontal tail root thickness ratio 

TC tt _ - vertical tail root thickness ratio 

CRCL, “ horizontal tail root chord, ft. 

HT 

CRCI^ - vertical tail root chord, ft. 

F and F ■ tail load parameters 

H V 

W_ * design gross weight, lbs. 

G 

B ht ■ horizontal tail span, ft. 

By.j, ■ vertical tail span, ft. 

B *» wing span, ft. 

Ely • fuselage length,, ft. 

SLM h - horizontal tail taper ratio 
SLMy ■ vertical tail taper ratio 

SKy ■ horizontal tail weight trend factor (See Figure v.1.9) 

SK TL " horizontal tail weight penalty factor (nominally » 1.0) 

SK ■ vertical tail weight trend factor (See Figure V.1.10) 
z 
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The trends consider the tail loads which are a function of the gross 
weight, span, radii of gyration, and point of load application. The 
"S " tern in the vertical tail equation accounts for T- tail configur- 
ations. Figures V.1.9 and V.1.10 present the aircraft used to develop 
the trends. The term "sk tl " is in namelist XNGASP and has a default 
value of 1.0. It is included to provide a means for penalizing the weight 
of the horizontal tail when additional design loads, such as carrier 
landings, are a consideration. SK^ would then have a value between 1.1 
and 1.2 depending on the magnitude of the design loads. 

The terms "SK 11 and "SK '* are the weight trend factors for the 
y z 

horizontal and vertical tails respectively, as defined in reference 1. 
These items are related to the pitch radius of gyration for the 
horizontal tail and the yaw radius of gyration for the vertical tail. 
Figures V.1.11 and V.1.12 show the definition of these terms in relation 
to the radii of gyration as shown in the VASCOMP program, Reference 1. 

This relationship is based on the regression analysis presented in 
Reference 2. Since this analysis mainly contained a luge airplane data 
base some smaller airplanes were analyzed using these weight trend 
equations. Figure V.1.13 shows the resulting SK^ and SK^ terms for these 
smaller airplanes. 

V. 1.2. 2. 7 Wing Weight 

The wing structural weight trend equation is based on a semi- 
empirical relationship developed from the approech outlined in Reference 
3 and a regression analysis of the 18 aircraft shown in Figure V.1.14. 
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VERTICAL TAIL WEIGHT - POUNDS 


9,000 
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PITCH RADIUS OF GYRATION 


100 


60 


60 

SO 

40 


30 


1. 707, DC-6B, 

2. B-26, B-29, 


DC-3, 

B-47, 


404, 

B-52, 

"‘ T "" 


ETC., SKy - 0.17 
ETC., SRy - 0.16 



20 


* 10 
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4 
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10 


20 
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FIGURE V.1.11 - RADIUS OF GYRATION, PITCH. 
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YAW RADIUS OF GYRATION 





Aircraft 

SK v 

SK 

Z 

Cessna 150 

0.204 

0.227 

Bede 5J 

0.14 

0.190 

Arrow 

0.0747 

0.22c 

Cessna 210 

0.115 

0.164 

Cessna 340 

0.063 

0.10 

Citation 

0.102 

0.200 

Lear jet 

0.180 

0.220 

DHC-6 

0.110 

0.220 

Gulf stream I 

0.166 

0.48 


FIGURE V.1.13 FACTORS COMPUTED FOR SMALLER AIRCRAFT 
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Aircraft 

Gross 

Weight 

lbs. 

Wing 

Area 

ft. 

Actual Total 
Wing Weight 
lbs . 

Predicted Total 
Wing Weight 
lbs. 

Actual Weight 
of High Lift 
Devices** li is. 

1. 

Bede Si 

1000 

38 

90 


4.2 

2. 

Cessna 150* 

1600 

157 

206 

199.88 

lo.« 

3. 

Cessna 172* 

2300 

174 

235 

246.34 

17.4 

4. 

Piper Arrow 

2600 

160 

273 

275.70 

15.(> 

5. 

Cessna 182* 

2800 

174 

257 

281.98 

16. < 

6. 

Cessna 210 

3800 

175 

375 

395.87 

22. i 

7. 

Cessna 340 

5975 

185 

544 

530.59 

17.2 

8. 

CITATION 

11650 

268 

1019 

1064.61 

85.(i 

9. 

DHC-6* 

12500 

420 

1212 

1480.17 

134. 1 

10. 

Merlin IV 

12500 

277 

1316 

1249.80 

100 

11. 

Lear jet 24 

13500 

232 

1175 

1151.22 

170 

12. 

Fregate 

23810 

600 

2749 

1981.87 

320 

13. 

Gulfstream I 

35100 

610 

3643 

4162.57 

472 

14. 

F-27 

35700 

754 

4265 

4121.93 

516 

15. 

Jetstar I 

40921 

543 

2954 

3914.50 

467 

16. 

F-28 

65000 

822 

7135 

7027.85 

928 

17. 

737-200 

100000 

1005 

10775 

10332.96 

2426 

18. 

DC-9- 32 

108000 

1001 

11391 

11409.55 

2241 


* STRUT- BRACED WINGS **includes trailing edge and loading edge devices 


FIGURE V.1.14 WING WEIGHT DATA 
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The method outlined in Reference 3 is based on a generalized 
expression for the material required to resist the root bending moment 
due to wing lift in a specified flight condition. The method does not 
include estimating the weight of the high-lift devices (see subroutine 
FLAPS) but does account for type of wing, wing relief factors, and 
landing gear location. The resulting trend equation is of the following 
form. 


10— S 1 044 4 

w w - SK *w sk no S1 W sk gear f oo b • ♦ S V ♦ w , 


where 


V" CS WC2 


1.535 


HLDEV 

(V.1.67) 


SK WW " w * n 9 weight trend factor in namelist INGASP 
(default ■ 133.4) 

SK nq - correction factor for the non-optimum material 
- 1 - 2 . 5 ^ ^T/ CS. 


(V.1.68) 


WC2 


(V.1.69) 

SK EP0S " en 9i n e position factor from Reference 3 

■ 1.0 propeller aircraft with no wing-mounted engines 

■ 0.98 propeller aircraft with 2 wing-mounted engines 

■ 0.95 propeller aircraft with 4 wing-mounted engines 

or high-subsonic jet aircraft with 2 wing-mounted engines 

■ 0.90 high-subsonic jet aircraft with 4 wing-mounted 
engines 

- 1.05 high-subsonic jet aircraft with no wing-mounted 
engines 


SK _ ■ landing gear location factor from Reference 3 
GEAR 

■ 1.0 for wing-mounted landing gear 

«■ 0.95 for landing gear not mounted on wing 


(V.1.70) 
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r oo • wln9 1 °* din » P““«“ 

- [ sk str u i* (“o - •" \ ) ] ™ ,V1 - 71 > 

SK - reduction in bending moment factor for Btreet braced 
5TR 

wing from Reference 3 

- 1 - (STRUT) 2 (V.1.72) 


STRUT • ratio of Bpanwise location of wing strut to wing 

semi -span (equals 0 for canti lever wing) 

U - ultimate design load factor, g's 

Lr 

W - design gross weight, lbs. 

G 

W ■ wing weight, lbs 

W 1 


SL 


M 


TC. 


CS 


WC2 


B - wing span, ft. 

- wing taper ratio 
* wing root thickness ratio 
-> cosine of wing half-chord sweeps angle 
W HLDEV " of wing high-lift devices, lbs. 

(computed in subroutine FLAPS) 

Figure V.1.14 compares the predicted and actual wing weights for the 
18 aircraft used in the regression analysis to determine the default value 
of SK^ and the exponents on span, taper ratio, thickness, and sweep 
angle terms. 


V . 1 . 2 . 2 . 8 Fuselage Weigh t 

The weight of the fuselage structure is determined from the weight 
trend equation developed in Reference 1 and presented below. 


w_ 


" sk b 


(k, 0 * 508 ♦ 


w, 


BOOM 


(V.1.73) 
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where 


k - (lO- 4 W x ) ■ (lO- 3 S p ) S w ( EW . EL „ ) ’ lo, (VD mis ) 

^DELp ♦ 1 ) 0 


.2 .3 

*LF 


(V.1.74) 


"BOOM “ 10 ‘ 3 SK BM EL BM ' /S * W W G 


(V.l. 75) 


w x 

S F 

S WF 

el ffc 

el rw 

VD 

MIN 

DEL.. 


LEGEND ; 

SK^ - fuselage weight trend factor In namelist INGASP 

(default ■ 136) 

* weight of fuselage and contents (includes empennage), lbs. 

2 

* wetted area of fuselage, ft. 

■ fuselage width, ft. 

■ length of fuselage, ft. 

■ length of pylon for fuselage mounted engines, ft. 

“ dive speed, kts 

- limit differential cabin pressure, psi 

■ ultimate load factor 

■ weight of fuselage ta:l boom, lbs. (KCONFG ■ 1) 

■ tail boom weight trend factor in namelist INGASP 

* length of tail boom, ft. 

■ cross-sectional area of tail boom, ft. 

■ aircraft design gross weight, lbs. 

v> 

The above trend equation was developed using a number df commercial, 
military, and cargo aircraft. Figure V.l. 15 indicates the relative body 
weight variation between these different families of aircraft. A mean 
line with an value of 124 is shown to represent the average of all 
the aircraft shown in Figure V.l. 15. Also shown is a body adjustment 
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BODY GROUP WEIGHT - POUNDS 
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1.4) MIL TYPE AIRCRAFT 
(MAX CUTOUTS, CARRIER USE) 

. i . : : i , i ■ ' 

1.)) NIL TYPE AIRCRAFT 

! i .1 I • I i i 

1.00 AVO SLOPE FOR ALL AIRCRAFT SHOW 

I ! ! I . I ?•:..! ; ’ 

[» » 0.04 AVC SLOPE TOR COWEPC1AL »lRCl(Arr I 

1 


Selection of the "a* feetor nay be detenlntd 
by examining the above aircraft and applying 
the a factor which beat describee the confly 
oration being sited. Conaideratlon oust be 
given to those Iters having a direct Influence 

r n the velght of a foselage, sueh as the rela- . 

iva number of cutouta, pressurlsatlon require-^ 
bents, reap site, etc. ror aircraft requiring 
additional structure over and above that which 
boreal ly exists in conventional designs (addi- 
tional fairings# supports, etet), the weight 
Density stay be accounted for In the *e* factor 
or in the incremental group weight or stultl- 
Dlieator factor sections of the weight input 
feheet. 
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FIGURE V. 1.15 - BODY GROUP WEIGHT TRENDS 
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factor "a" which corrects the 124 constant in accordance with the specific 
family being considered. The revised constant ( 124 M a n , is the SK^ term 
to be input in namelist INGASP. The limit differential cabin pressure 
(DELp) is input also in INGASP or the program calculates the appropriate 
value to provide an 8000 ft. cabin altitude at the design cruise altitude. 
Since the analysis of Reference 1 mainly contained a large airplane data 
base some smaller airplanes were analyzed using the weight trend equation 
above. Figure V.1.16 shows the resulting SKg for these smaller airplanes. 

V. 1.2. 2.9 Flight Controls Group Weights 

The weight of the total flight controls group is estimated first and 
then the weight of the cockpit controls component. The difference between 
these two weights is considered the fixed wing surface control portion of 
the flight controls. 

The weight of the total flight controls group is estimated from the 
following trend equation which is based on a number of military and 
commercial aircraft. 


“sex ■ s w - 317 <1°' 3 V 602 


(U )‘ 525 (Q ) 
' LF ' W DIV 


.345 


(V. 1 .76) 


where 

W ■ weight of total flight controls group, lbs. 

SKp w ■ flight controls weight trend factor in namelist INGASP 

(default ■ 0.404) 

2 

S - wing planform area, ft. 
w 

W - aircraft design gross weight, lbs. 

G 

U^p ■ design ultimate load factor, g’s 

2 

Q djv • dynamic pressure at design dive speed, lbs. /ft. 
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Aircraft 

Pressurized 

■ 

Cessna 150 

NO 

130.0 

Bede 5J 

NO 

160.0 

Arrow 

NO 

128.0 

Cessna 210 

NO 

161.5 

Cessna 340 

Yes 

101.0 

Citation 

Yes 

87.0 

Lear jet 

Yes 

136.0 

DHC-6 

No 

128.5 

Gulfstream I 

Yes 

128.5 


FIGURE V.1.16 COMPUTED FUSELAGE WEIGHT FACTORS 
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The cockpit controls weight is estimated as the following function 
of gross weight 


fcl - SK <l(f V )* 41 

CC CC G 


(V.1.77) 


where 


SK ■ cockpit controls weight trend factor in namelist XNGASP 
cc 


(default * 11.0) 

The fixed wing surface controls portion is then the difference 


wc - w 
FW SCX 



(V.1.78) 


Figure v.1.17 indicates seme typical values of SK^ and SK^ for a 
selected number of small airplanes. 

The sum of the aircraft flight control component weights is 


w 


FC 


CK, c Vi + CK, _ WC + CK, Vi „ + DELW 
15 CC 18 FW 19 SAS FC 


(V.1.79) 


where 

W - SK - weight of stability augmentation system input in 

bAb bAb 

namelist INGASP (default - 0.), lbs. 

DELW Fc ■ incremental flight controls weight, lbs., (this may 
be input in namelist INGASP) 

and where the constants CKj are currently the value of unity. If the 
flight controls weight trends are to be by-passed the user can input 
SIC, , SK , SK as zero agd input the total flight controls weight as 

r W LC bAb 

DELW . 

FC 

V. 1.2. 2. 10 Optional Structural Weight Computation Bypass 
If the structural weight trends are to be by-passed the user can 

input SK^, SK y« SK Z » SK 3 ' SK lg' SK PES 33 zero anc * i n P ut the total 
structural weight as DELW^. 
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Aircraft 

SK rw 

SK cc 

Cessna 150 

0.485 

11.5 

Bede 5J 

0.133 

3.0 

Arrovi 

0.383 

13.5 

Cessna 210 

0.322 

8.7 

Cessna 340 

0.350 

9.8 

Citation 

0.430 

20.6 

Lear jet 

0.404 

11.0 

DHC-6 

0.430 

20.0 

Gulf stream I 

0.430 

20.0 


FIGURE V.1.17 computed flight controls weight factors 
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where 


DELW gT - incremental etructural weight, lbs. (this may be input 
in namelist INGASP) 

V. 1.2. 2. 11 Fuel System and Fuel Weight 

The fuel system weight depends on the amount of fuel available (WFADES) 
for the design payload condition of the airplane. The design payload is 
input as WPLX in namelist INGASP and is identified as WPLDES in subroutine 
WGHT. If no value is input than WPLDES is a function of the number of 
passengers (PAX) and their unit weight (UHPAX) . 

The weight of the fuel system is basically a percentage of the fuel 
volume at the design payload condition as expressed by 

VJ „ - CK_, (6.687/FUEL ) SK (V.l.So) 

FSS 21 D FS FADES 

Since the design fuel load is determined by subtracting the component 

weights and payload from the gross weight, the design fuel load is 

expressed as 


w, 


FADES 


(W„ - W - W_ - - W„,„ - W „)/D (V.1.81) 

G PSTAR ST FC FE FUL PLDES 


where 

D * 1 + CK 21 SK PS 6 - 687/7UE S (V. 1.82 

W FUL " ^ xe< * useful load weight 

FUELp * fuel density in lb. per gal. (program sets it to 6.687 for 
jet fuel and 6.0 for gasoline) 

CK 2i *■ weight multiplication factor for fuel system, currently sot 
to 1.0 in program 

SKpg « .0195, default value of weight coefficient for fuel system 
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The value of SKy g will vary dapandlng on the capacity, type, and 

complexity of tha aystam required. For commercial aircraft having simple 

fuel systems in tha wing, the value for SK would range between 0.02 and 

FS 

0.07} for aircraft requiring self-sealing tanks with more complex systems, 
the value would range between 0.10 and 0.1S. Figure V.1.18 contains the 
value of SKpg for some general aviation aircraft. 

Tha total propulsion weight is than computed as 


w „ ■ ♦ W_„ 

P PSTAR FSS 


(V.1.83) 


where the first term har n defined in Bq. (V.1.55) and the second is 
the weight of the fue a\ ^xy system from equation V.1.80. 

In addition to the design fuel load, two other fuel loading conditions 
are calculated for determining the payload-range calculations performed 
in XRANGE. One is the minimum fuel condition (WFAMIN) and the other 
the maximum fuel case (WFAMAX) . 


M gg U 0 U mm U M U a U a U w U 

FAMIN G PSTAR ST FC FE FUL FSS 


- W, 


PLMAX 


(V. 1 .84) 


The W FAMAX equal to either w FADE r or tlie * uel available for filling 
all the available volume without exceeding the gross weight. 

The next group of statements deal with determining the available wing 
fuel volume and weight. The wing tank fuel volume is found in cu. ft. as 


VOLW 


.8889 SK^ TC S*' 5 (2.*SL M ♦ 1) 

V^R (SL + i ) 2 
n 


(V.1.B5) 


where TC is the average wing thickness to chord ratio, calculated on the 
previous line. SIC— is a factor that relates wing volume available for 
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Aircraft 

sk fs 

Bede 5J 

.050 

CE 150 

.119 

Piper Arrow 

.138 

CE 210 

.0394 

CE 340 

.0415 

CE Citation 

.052 

Twin Otter 

.060 

Lear jet 

.020 

Gulfstream I 

.0143 


FIGURE V.1.18 CALCULATED FUEL SYSTEM WEIGHT COEFFICIENTS 


V-l 

47 






(ual to total theoretical wing volume, xf SK^ la input as aero than no 
fuel will be carried in the wing. For completely wet winge SK^ it on 
the order of .43 to .45. Figure V.1.19 shows the value of SK^p for some 
general aviation aircraft. The default value of SK^, is .43. 

The logic between statements 38 and 60 sets up an iteration to 
converge on a wing weight since fuel weight can be considered as a relief 
load factor. 

V.l.2.3 Tip Tank Siring and Weight 

After the wing weight iteration has converged* and if the maximum 
available fuel is more than the wing fuel capacity \ W__) * the 

tip-tank sizing calculations will start if that option has been selected 
(KTIP * 1) . The tip-tank calculations determine the size of two tip 
tanks containing the oxcess fuel. 

The tip tanks are simulated as prolate spheriods having a major 

axis* A , and minor axis, B . The dimensions of these axes must 
XXS XI5 

satisfy the constraints shown in V.1.20. The tip-tank sizing starts by 
sizing a tank based on the minimum tank length and growing the length by 
increments of .25 ft. until all the extra fuel is accomodated or the 
maximum allowable tank size is reached. Constraints are put on tank 
length, diameter, and fineness ratio (L/D) . The si>: constraints are 
illustrated in Figure V.1.21, in which the length, diameter and volume 
vary as indicated. For > L i / D 2 and R 2 ( L 2 /,D l' * as shown * * two o£ thu 
extreme corners of the rectangle cannot be used. Notice that the locus of 
fixed tank volume may not contact any of the four limiting values of 
length or dit-.ter (locus V*). 
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Aircraft 

S1 W 

Remarks 

Bede 5J 

0.43 

wet wing 

CE 150 

0.0762 

only uses portion of available volume 

Arrow 

0.0783 

only uses portion of available volume 

CE 340 

0.114 

majority of fuel carried in tip tanks 

Twin Otter 

0 . 

Fuselage fuel tanka 

Lear jet 

0.43 

wet wing plus tip tanks 

Gulf stream I 

0.3822 

wet wing 


FIGURE V.1.19 CALCULATED WING FUEL VOLUME COEFFICIENTS 
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■n B xis I" - 





The tip tank volume in cu. ft. ie determined from the tip tank fuel 
requirement end the fuel density, FUEL^. 


V mD - W / (7.4805 * FUEL ) (V.1.87) 

FTP FTP D 

The wing root chord is derived from the mean chord (CBAR^) and the taper 
ratio (SL m > as 

CH r - | CflAiy' [l ♦ SL m - Sly' (1 ♦ SI^)] (V.1.88) 

and the tip chord ie 


CH T " SL M CH R (V.1.89) 
Limiting values for tip tank dimensions are defined, in terms of which the 
minimum and maximum values of tank volume are given by 


V FTPMN ■ S Vr “MIN “MIN 

and 2 

(1.046) AX BX^ 
FTPMX FT MAX MAX 


(V.1.90) 

(V.1.91) 


In these expressions for volume of a prolate spheroid, the numerical 
factor is 2. (*/6.) to account for both tanks, and the factor SK^ is 
input at .979 to account for volume of the tank that is not structure. 

The tank diameter needed for minimum length is found, if the volume 
exceeds the minimum volume, as by inverting the volume equation; i.e., 


B XIS ■ t' 956 V FTP / ‘“is • SI W ] 1/2 < V - 1 - 92 ’ 

If this is within bounds, control is transferred to statement 30, where 
succeeding calculations include an equation for total surface area of the 
tip tanks. For two prolate spheroid tanks, the total area is found as 
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S TIP “ ® ^ B xis * *xis B xis 8in ^XSS^XIS^ 

(V.1.93) 

’ 2 * 11 * ^IN "W 4 * 


where 

E 


XIS 


-Cl- 


<B X1S /A X1S ) 


2 T 1/2 


(V.1.94) 


E Is the eccentricity of the ellipse of revolution. The last term in 
xxs 

equation V.1.93 approximates the area of the wing tip as that of an 
ellipse with axes &X MIN and bx mjn' and is subtracted because it is not 
part of the wetted area of the tip tank. 

The empty weight of the tanks is proportional to the surface area 


“tip ■ SI W 5 tip (v - 1 - 95 > 
where SK^p - 1.89 is input as the weight in lb. per sq. ft. of the tank. 
The total weight of tip tank structure and fuel is 


W C " W TIP + W FTP 


(V.1.96) 


and the remaining equations in this portion of the program are concerned 
with the iterative tip tank parameters contained between statements 
numbered 5 and 80. 

V.1.3 Aircraft Balance 

During the task of concept formulation and aircraft definition, 
aircraft balance can be a major concern. Th*s problem does not stand by 
itself. It reaches into other aspects of aircraft design such as stability, 
control, performance, structure weight, engine location, payload require- 
ments, and cost. Two balance options are currently provided for in GASP. 
The first (LOWING - 1) is a simplified method and the uecond (LOWING = 2 ) 
goes into more detail stability, control, and loading requirements. 
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V. 1.3.1 Simplified Balance Option 

The simplified method is based on the following conditions i 

(a) Horizontal tail volume coefficient is specified and held 
constant; i.e., the horizontal tail area varies inversely 
to the tail length. 

(b) First-order values for wing-plus-body aerodynamic center, 
horizontal tail contribution, and static margin must be 
known. 

(c) The airplane shall be balanced with its center of gravity on 
a certain point on the wing mean aerodynamic chord. 

Following the tip tank calculations, the program deals with the 
relocation of the wing if the balance option flag is LC WING ■ 1 or 2. 
This portion begins with the computation of remaining weight, the 

sum of fuselage, payload, fixed equipment and controls weights. This is 
found by subtracting many component weights from the gross weight, or 


rj B M • W • u . u - W "*W - w - w 

REMN G HT VT W PSTAR FW FTP TIP 


(V .1.97) 


- - CK_ V* - XW _ 

PES BOOM 21 FSS LG 


This is followed by a sequence of geometrically derivable quanti- 
ties, related to wing, horizontal tail and vertical tail. For the wing, 
these include the wing centerline chord. 


CR CLW* 2 - VC E (1 + S Vl 
the tangent of the leading edge sweep angle. 


T 


SWPLE 


TAN (RLM c4 ) + (1 - SL m )/ [AR (1 + SL^)] 


(V.1.98) 


(V.I.99) 
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the mean aerodynamic chord. 


*MAC ' 8 * T S«PLE <l + 2SL H ,/ [*• (1 * SL M>] (V. 1.100) 

and the distance from the aircraft nose to the wing centerline quarter chord, 

EL C4W " EL F ' CR CLVT + ’"'mAC ’ EL TV “ Vc 

(V. 1.101) 

+ .25 [CBAR^ - CBAI^ 4 CR^J 

where SL is the wing taper ratio, and where the various distances in 

M 

V.I.101 are shown in Figure V.1.22. The root chord and mean aerodynamic 
chord are also developed for the vertical and horizontal tail surfaces, 
using analcgous equations and symbols. Following these are expressions 
for the tip chord, which duplicates the equation found in Eq. (V.1.89) 

CH T “ | SL m CBAR^ [l + SL M - SL m / (1 + SL m >] (V. 1.102) 

and for the distance from the aircraft nose to the half chord of the 
wingtip. This latter distance is given the name 

EL " EL - . 2SCR + x _ + . 25CBAR (V. 1.103) 

TIP C4W CLW MAC W 

If the wing tip tanks are centered here, their trailing edge occurs at the 

distance 


EL - EL + .5 A , 

TIPE TIP XIS' 

aft of the nose. The wing tip trailing edge is located at 

EL - EL + .5 * B * TAN (RLM ) + .75 C , 
TPMX C4W C4 HT 


(V. 1.104) 


(V. 1.105) 


and if EL < EL m n «v' aB occurs for a swept-back wing, the centers of 
TIPE TPMX 

the tip tanks are located 1/2 tank length ahead of the tip trailing edge. 
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The following computations deal with the moment arms from aircraft 
nose to the estimated center of gravity of vertical and horizontal tails, 

EL COV ■ EL C4W - -““CLW * Vc * - 25CaW> W * E Sv * • 08CBA V 
el cgh ■ el cgv - w •“ (cba *vt - '“Sr 1 ♦ W 

♦ .25 (CBAK^J, - CKCL^.) ♦ El I(|c ♦ S \ V < S “ P fi C v> (V.X.107) 

where the last term in Eq. V. 1.106 assumes the c.g. of the vertical tail is 
at 8« aft of the mac 1/4 chord point for vertical tail. The other undefined 
quantities in this equation are 

■ vertical location of horizontal tail on vertical tail, 

0 < SA„ < i. 

H 

■ span of vertical tail 

- aft displacement of horizontal tail root chord leading edge 
relative to vertical tail root chord leading edge 
The program follows with equations for the spanwise locations of the 
mean aerodynamic chord of the wing, 

Y MAC “ B (1 + 2S V/ t 6 (1 + s yJ (V. 1.108) 

and a similar expression for the horizontal tail mac location, where SL 

M 

is the taper ratio of the wing. 

For the simplified option “ D the program now passes to 

statement 100 where the aircraft center of gravity location aft of the 
nose is then expressed as. 


sa h 


b vt 

EL 


INC 
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el cg ■ ['Shop <m psta* * w pes ) * CBAR w (cp mrgn ~ st krgn ) 

<"« *%*»«* ck 21 m fss> * 'Sip ‘"tip * W 

(V • X • lv7 J 

* el cgh w ht + el cgv w vt * el remtj m bemm * "book 

< E S - • “Sm 1 ]/ < W o - S - M p« - ™V 3 - <*21 "pss> 


Zn terms of this distance, the wing center of gravity is the follow- 
ing distance aft of the nose. 


E W - El C0 ♦ (CP MM» • S W 


(V. 1.110) 


where 


CP un _, ■ .10, location of wing CG in relation to the quarter cord 

MRuN 

of the wing mac 

ST„_, ■ 0., location of aircraft CG in relation to the quarter 

MRGN 

cord of the wing mac 

The distance from nose to the quarter chord of the wing mac is 


EL CP “ el cg - st mrgn cbar W 


(V. 1.111) 


This equation (Statement 102) is also the point in the program where 
control returns after the more detail balance option. The next dozen 
equations are related to parameters which are descriptive of the horizontal 
and vertical tail. These include the moment arms for both tail surfaces, 
and these equations can be derived from study of Figure V.1.23. 


el tv “ el f " el wing ” cr clvt 


♦ ♦ * 25 CBAR VT 


+ CP 


MRGN 


CBAR W 


(V. 1.112) 
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25 (CBAS^ ♦ CRCI^ - CBAA^) 


(V. 1.113) 


EL - EL 
TH TV 


+ XH . 4- EL ♦ SA B tan (SWP ) 

MAC INC H VT QCV' 

The areas of the horizontal and vertical tail surfaces are then re- 
calculated 

S ht “ VBAR H s w CBA V e Sh (V. 1.114) 

S VT - ^ S W B/E Sv (V.1.U5) 


and VBAR^ are dimensionless tail volume co- 
efficients. Other geometric calculations are concerned with span, root 
chord and mean aerodynamic chord of both tail surfaces; using straight- 
forward geometric equations which have already been discussed. 

The criterion for convergence in the tail sizing is that the dif- 
ference in successive horizontal tail areas be less than .4 sq.ft. 
Otherwise, control returns to statement 4, where the tail sizing loop 
begins. 

At statement 108, the cruise altitude, Mach number and weight are 
redefined, and subroutine CTAER is called to get the drag of the resized 
configuration. The aircraft center of gravity locaion on the wing mac is 


where the parameters VBAR^ 


°CGMAC " 0 L CO * EL CP * CBAR «] / CBA *W (V.1.U6 

Following a number of written statements regarding geometrical parameters, 
the subroutine then computes: 

(a) leading edge sweep angles of wing, horizontal tail and vertical 

tail, in degrees, SWP_ . SWH , SWV , 

LED LED LED 
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(b) tip chorda of horizontal and vortical tails CH^# CH^, 
after which 28 more geometrical parameters ara printed. 

The operating weight empty is determined as the svai of the weights 
of the propulsion group, flight controls, structures, fixed equipment and 
fixed useful load is then given by 


O t „ - W ♦ v ♦ u ♦ v * v 
WE P FC ST PE FUL 


(V. 1.117) 


and if the sum of this and the maximum available fuel weight exceeds the 


gross weight W Q , the maximum available fuel weight is redefined as 

W a U ~ O 

FAMAX G WE 


(V. 1.118) 


Now, if I ) 0, or ■ 0 or - 2, the program returns. 

These flags correspond to no tip tanks being added or engine resized, in 
which case the drag and engine size do not require recomputation. Other- 
wise, the subroutines CTAER, ENGSZ and ENGWT are called prior to the 
second iteration of the loop beginning with the minimum wing fuel weight. 

V.l.3.2 Detail Balance Option 

The detail balance option tries to position the wing so that the c.g. 
travel from loading considerations falling within the allowable c.g. 
range from stability and control considerations. The logic for this 
option is as follows: 

(1) The wing is initially positioned based on input or default 
tail volume coefficients and forward and aft stability limits 
are assumed . 

(2) Based on c.g. limits the most forward and aft c.g. positions 
from a loading standpoint as well as the c.g. travel are 
determined. The loading rules are summarized later. 
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(3) Rs-position wing so that c.g. at most forward loading condition 
(ELCGF) falls at forward c.g. stability and control lljait 
(FCGLIM) . Compute tail moment arms . 

(4) Call subroutine TAIL with c.g. travel required to compute new 
tall areas and volumes and new stability and control limits. 

(5) Compare tail sise as in simplified balance option, if not 
within allowable tolerance repeat process starting at (1) 
above. 

The group of statements in WGHT starting at the comment "COMPUTE 
MOST FWD AND AFT LOAD CONDITIONS" through statement 100 deal with the 
computation of the extreme center of gravity locations due to most for- 
ward and aft load conditions. These conditions are based on the loading 
rules summarized in the following sections V.l.3.3 to V. 1.3. 3. 7. 

FCGLIM refers to the forward c.g. limit and ACGLXM is the aft c.g. 
limit from stability and control considerations. 

The forward and aft loading conditions are specified in terms of 
the extreme values of center of gravity, expressed as a percent of mac. 

The forward c.g. location, in ft., is 

F CGLIM “ EL WING ’ (CP MRGN + * 25 ” *CGWD^ CBA *W (V. 1.119) 

and the aft location A__. is of identical form, where the parameters 

CGLIM 

X CGFWD and x cgaft are 9 iven as percent mac. The center of gravity for 
the operating weight empty condition is 

CG^ itc , - EMO /O (V. 1.120) 

OWE WE WE 

where the two factors are defined in terms of specified weights and 
moment arms for major aircraft components. The number of passengers, 

N , and the number of rows of seats, N , are then found under the 

PAX KUW 
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assumption that each passenger weighs tha input unit weight, with baggage, 
in tanna of tha number of aaata abreast, S^. 

The loop 81 ia concerned with tha moment produced by rowa of passen- 
ger a, measured relative to the nose of the aircraft. Thus, the moment 
arm is, for row nutobir I, 

“•KW «> - V l *-‘ 1 * el odn m n * el pc ■ V 36 - (V.l.ui) 

where the last three terms measure from the aircraft nose to the front of 
the fr^raenger compartment » i.e., 

E Lqdn “ * in ® n ® 88 ratio of fuselage nose 

- height of fuselage nose, ft. 

EL pc - length of pilot's compartment, ft. 

P ■ lengthwise distance between sects, in. 

5 

This permits the computation of passenger-induced moment as the sum, 

EM "EM ♦ s (170) EL (I) (V. 1.122) 

PAXT PAXT AB ' ’ ROW 

The more complex loop 88 is intended to develop the moment caused 

by several component weights when the c.g. is at the most aft (X * 1), 

the most forward ■ 2) and the design location (X. TU __ • 3). The 

LIMIT LIMIT 

loop is called three times, corresponding to these c.g. positions, and 
the computations are complex chiefly because of the number of component 
weights which can vary on the aircraft. 

The first third of the loop is devoted to computation of total weight 
and moment of the passengers, the moment being measured relative to the 

nose of the aircraft. These are denoted by W PAXD ^ LIMIT ^ anc * ^paxd 

(I ) , respectively, and are found by straightforward means, being 
LIMIT 
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proportional to number of seats abreast, S__ and to unit weight of 

AB 

passengers, UW pAJ{ . 

Following this segment, the baggage weight and moment arm are found 

“ "bag "limit 1 “ 4 el bags ^ o£ p ° ylold ' P“ sen9er “ l « ht and 
proportion of fuselage length, ELp. The fuel weight airt center of 

gravity are then found in terms of wing geometry and kito^n limits of 

center of gravity location. Tip tanks are included in this computation, 

which ends w ' v h an expression for 


X_ IIE<r <I T (V. 1.123) 

TUELD LIMIT FUELM FTPD LIMIT FWD LIMIT 


where I “1 or 2 • and 

LIMIT ' 


“fUELD (3) “ W 3 “ °WE - W PAXD (3 > “ W BAG (3) * W FWD < 3) 


(V. 1.124) 


- W____ (3) 
FTPD 


These weights must be positive, and they contribute to fuel weight, 

W FUEL “ W FWD ^ LIMIT ^ + W FTPD ^LIMIT* + ^FUELD ^ LIMIT ^ (V. 1.125) 

and hence to the limit weight, 


W T ^LIMIT* “ °WE + W FUEL + W PLD 


(V. 1.126) 


where 


WE 


■ operating empty aircraft weight; lb. 


W CTtT> ■ wing fuel weight, lb. 

FWD LIMIT 

WFTPO (1 LIMIT > " tiP tank fUel Wei * hfc ' lb * 


W 


PLD 


■ weight of pay load, lb. 


If the limit weight exceeds W_, baggage weight is redefined as 

G 


v-i 

64 



W BAG ^LIMJT* " W BAG (I LIMIT ) ” W T (I LIMIT ) + W G 


(V. 1.127) 


and 


W T ^LIMIT* “ W G 


Flnally« the center of gravity ia located the following distance aft, 


EL (I ) ■ EM (I )/W (1 ) 

CGD ' LIMIT' L LIMIT" T ' LIMIT 


(V. 1.128) 


where the divisor is in Eq. V. 1.126, and where the moment ia the sum of 
five terms. 


^L (I LIMIT } “ ^OWE 4 EL SING [ W FWD ^ LIMIT ^ + *FUELD 

<I LIMIT > ^ + EM PAXD (I LIMIT > + EL BAGS W BAG (V. 1.129) 

^LIMIT* 4 E SlP W FTPD ^LIMIT^ 


These terms correspond respectively to operating empty weight, wing fuel 
weight, passenger weight, baggage weight and tip tank fuel weight, for 
the three loading conditions, I TTUT _ * 1, 2 or 3. 

LIMIT 

A large number of parameters are printed, following which c.g. moment 
arms and wing 1/4 chord locations are calculated as EL cp* EL CG and 

G CGMAC' in termS ° f CP MRGN' X CGFWD' CBA *W and related terms * 

A call is next made to .subroutine TAIL, which is described in 
section v.1.4, to determine tail size for the required c.g. travel, 

DELXCG . 

V.l.3.3 Most Forward Load Condition 

1. Pilot 

2. All seats forward of FCGLIM are £1 , by passengers each 
weighing 170 lbs. 
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3. If baggage lies forward of FCGLIM, all baggage la loaded, 
otherwise none (baggage load defined below) 

4. All fuel tanks forward of FCGLIM s full 

5. If fuel forward of FCGLIM is less than minimum fuel (see below) , 
add fuel to next most forward fuel tank until minimum fuel load 
is reached 

V. 1.3.4 Most Aft Load Condition 

1. Pilot 

2. All seats aft of ACGLIM are filled by passengers each weighing 
170 lbs. 

3. If baggage lies aft of ACGLIM, all baggage is loaded, otherwise 
none 

4. All fuel tanks aft of ACGLIM are full 

5. If fuel aft of ACGLIM is less than minimum fuel, add fuel to 
next most aft fuel tanks until minimum fuel load is reached 

V.l.3.5 Design Load Condition 

1. Pilot 

2. Number of passengers is largest whole number contained in 
WPL/200. Passengers are loaded in forward seats 

3. Baggage is WPL minus number of passengers x 170 

4. Design fuel load contained in wing, tip tanks (if any) , and 
fuselage tanks (if any) 

V.l.3.6 Minimum Fuel for Weight and Balance* Condition 

1. Piston engine aircraft: Minimum fuel (lbs) = maximum installed horsepower 

2 

*Min fuel for wt. and balance is not necessarily the same as minimum fuel 
for flight. See FAA Advisory Circular AC65-9, p. 50-51. 
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2. Turbine Aircraft! Specified by Manufacturer (no general rule) 

The program currently assumes 20% of wing 
fuel tank capacity 

V.l.3.7 Maximum Baggage Condition 

1. Maximum baggage is assumed to be same as design load baggage 
(see above) . Baggage is located by input variable RELB 

V.1.4 Stability and Control (TAIL) 

Subroutine TAIL is a computerized methc/dology used to determine the 
critical longitudinal and lateral stability and control requirements 
governing horizontal and vertical tail sizes of aircraft from two place 
trainers to business jets. 

The procedures which were developed reflect the requirements of the 
Federal Aviation Regulations (FAR) and accepted engineering design 
practices. The tail sizing procedures account for the following: 
Horizontal 

(a) Meeting a specified static margin in relation to neutral static 
stability (dcm/dc ■ 0) 

L 

(b) The ability to rotate the aircraft to takeoff 
rotation speed 

(c) The ability to obtain the maximum lift coefficient of the 
aircraft in the landing configuration in grovnd effect. 

Vertical 

(a) Static directional stability requirements (weathercock 
stability) 

(b) Minimum control speed requirements for twin-engine aricraft. 
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To best understand the computations in this subroutine is to first 
review the output from the subroutine which summarizes the Important 
tail sizing parameters. Figure V.1.24 is an example of this output. 

The first three lines of print are the parameters that depend on *he 
particular flight condition for horizontal tail sizing. These are 
cruise, nose-wheel liftoff, and landing flare. The parameters are wing 
angle of attack (ALPHA) , wing lift curve slope (WING CLA) , tail lift 
curve slope (TAIL CLA) , tail efficiency factor (TAIL EFF) , downwash 
angle at tail (DOWNWASH) , wing lift coefficient (WING CL) , and moment 
derivatives (DCM) or moments (CM) for the fuselage, nacelle, flap and 
power effects. The next four lines of print deal with the elevator 
characteristics and effectiveness as identified on the output. The next 
five lines display the critical c.g. limits and the horizontal tail sizes 
required for the different criterion. The neutral point, aft and forward 
c.g. limits, and various tail sizes are determined in TAIL. The other 
items have been input or determined in WGHT. The last three lines of 
print are the vertical tail sizing requirements and the required vertical 
tail area. 

Initially the subroutine contains about 60 COMMON and DATA statements 
that transfer data and store data for table interpolations. The last 
data statement is important in that it makes the assumption that the 
aerodynamic center of the wing with flaps up (XAC) and flaps extended 
(XACF) is at the quarter chord. The next statement sets up a function 
statement for calculating the lift curve slope of the wing or tail 
knowing the aspect ratio, sweep, and Mach number. This is the same 
method that is used in subroutine CLIFT. 
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FIGURE V.1.24 - TAIL SIZING OUTPUT 


V. 1.4.1 Horizontal Tall Sizing 

After some initialization statements a number of parameters are 
determined for the cruise condition. These are the cruise lift 
coefficient (CLCR) , cruise angle of attack (ALFC) and lift curve slope 
(CLALCR) from CLIFT, and for propeller aircraft the thrust coefficient. 


THRUST 


'TCR 


2 2 

1.4 EM D* „ 
0 PROP 


(V. 1.131) 


where P Q is static pressure and D pRop is propeller diameter. Next, is a 
cruise value for tail efficiency, 


^CR " [* 9 + * 1S Ah 3 t 1 * * 8 * C TRC /3 * 14159 3 (V. 1.132) 

where * position of horizontal tail on vertical tail. 

The next portion of TAIL deals with conditions at nose-wheel liftoff 
for a speed 05 percent of the stall speed in takeoff configuration. For 
propeller aircraft this includes the thrust coefficient and tail 
efficiency calculations as in the preceding cruise portion. 

The next section determines the "minimum control speed" and correspond- 
ing thrust used in the vertical tail sizing of multi-engined aircraft. 

The "minimum control speed" is an input fraction of the stall speed in 
the takeoff configuration; * 


MC 


FV V 
MCS SOGE 


(V. 1.133) 


The corresponding thrust, as returned by ENGINE, is used in calcu 
lating the thrust with one engine failed 


F NMC “ THRUST (EN P " 1 -> 


(V. 1.134) 
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and for propeller aircraft# this is used in defining the thrust coefficient, 
C TMC * And tail efficiency# ETA^, on the next lines. 

If the configuration utilizes turbo jet/turbofan engines# the next 
section computes the tail efficiency factors for the four flight conditions 
as a function of the horizontal tail position# SAH. 

The position of the main gear from the aircraft center of gravity 
is determined 

X G “ (CX A “ GCG mAC ) CBAR W (V. 1.135) 

where - distance of main gear from aircraft c.g.« ft. 

G 

CX - location of gear on mac, fraction of mac 

GCG - aircraft c.g. location on mac# fraction of mac 
MAC 

The following statements deal with determining stability factors for the 
fuselage# Q FUS » and the nacelle, Q NAC » or, Q N depending on whether the 
nacelle is on the wing or fuselage. The nacelle factors are only 
determined for multi-engine aircraft. 

For the fuselage and wing mounted nacelles, Q FUS and Q NAC are determined 
from Figure V.1.25 where X^ R is the position of the wing root 1/4 chord on 
the fuselage and X 4RN is the position of the wing root 1/4 chord on the 
nacelle. These are determined by 

X 4R‘ [ EL C4W* “V 2 - 1 ™ < du) C 4>] /EL F (V. 1.136) 

X 4RH ‘ On ' - 7S C V C - TSW te “"V 2 -! /EL » (V.1.137) 

" here CR NAC ' C RCLW C 1 - - < X - - S '„> "W <B/2 'd (V.1.130) 


V-l 

71 



FIGURE V.1.25 - FUSELAGE/NACELLE STABILITY FACTOR (FROM NACA TR711) 



■•call* naaiaua dlutttr 

FIGURE V.1.26 - NACELLE FINENESS-RATIO CORRECTION FACTOR 

1 
2 





Yp B/2 - DBAR^j/2 


(V.1.39) 


YN 


ACR 


C D _ T „ • wing chord at wing canterline, ft. 

RCLW 


SL ■ wing taper ratio 
M 


m fuselage width, ft. 


EL ■ nacelle length, ft. 
N 




nacelle diameter, ft. 


TSW - the tangent of the wing trailing edge sweep 
TE 

If the nacelles are located on the fuselage the stability factor. 


Q n , is determined from Figure V.1.26 as a function of nacelle fineness 
ratio. 


The next parameter determined in TAIL is the downwash gradient at 
the horizontal tail. This is identified as OEDAW and is based on an 
empirical method of estimating the downwash gradient behind straight- 
tapered wings given in DATCOM. 


DEDAW 


CL 


AEMCRU 


CL 


AO 


4,44 C r AR F SLM F HT P LAM 


1.19 


(V. 1.140) 


The factors F„_, F , F, . P f , u are wing aspect ratio, wing taper ratio, 
AR SijM HT LAM 

horizontal- tail-location, and wing sweep factors, respectively, determined 
from 


F - l./AR - 1./ (1. + AR 1 * 7 ) 

aK 


"su. ■ (1 °- - 3 s lm ,/7 - 


f ht ■ u - ■ V B,/ <2 el th / b> 

1/2 

F “ COS (DLM „ ) 

LAM ' C4 


1/3 


(V. 1.141) 
(V. 1.142) 
(V. 1.143) 


COS (DLM) 
C4 
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(V. 1.144) 



where AR, S^, DZiN C 4 and B u * wing aspect ratio, taper ratio, 1/4 chord 
sweep, and span, respectively, is the height of the tail relative to 
the wing chord plane and EL^ is the distance between the 1/4 chord of 
the wing mac. and the 1/4 chord of the tail mac. The ratio of the wing 
lift- curve slopes of cruise Mach number and zero Mach number, 

CL AEMCRl/ cl ao ) a PP roximates the M 34 * number effect on the downwash. 


Fuselage and nacelle stability contributions, dC /dc , are determined 

ED L 

in the next few program statements. They are functions of geometry and 
the factors Qp US » $ NAC and C N » which were previously determined. DCM p us 
and are the fuselage and nacelle contributions at cruise. 

DCMpst and DQIj^ct are the fuselage and nacelle contributions at nose- 
wheel liftoff. DCMp SL ^ are the re8 P ective contributions at 

landing conditions. The fuselage and nacelle pitching moments are then 
computed for the liftoff and landing condition. O^uqt and CMp rj being 
the fuselage pitching moments and CM NACT and CM^ t * ie nace ^^ e Pitching 
moments . 

The next items of importance are the hinge-moment coefficients and 

the elevator power and effectiveness. The two hinge-moment coefficients 

are, C , the "floating tendency", and, C ,, the "restoring 
hALF j HDEL3 

tendency". They are computed in terms of section parameters C Harp and 

C from Figure V.1.27 in the following equations 
HDEL 



C HALF3 

• c c /c 

HALF LALFT' LALTO 

(V. 1.145) 


C HDEL3 

“ C HDEL + TAU H (C HALF3 “ C HALF > 

(V. 1.146) 

where 

C LALFT 

- lift curve slope of horizontal tail 
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FIGURE V.1.27. 


. Section hinge moment parameter* vermt flap chord /airfoil chord 
ratto (\A(‘.V Onoi» an Fmm V \C V Wit L-ftil, '‘Wind tunnel Data on tha 

Arrwij nani.c ClunfimtliM of VirpUne Control Surface*," by R. 1. Seva. 



C. ■ lift curve elope of horizontal tell with 
LALTO 

infinite aepect ratio 

TAU ■ elevator effectiveness 
H 

The next few lines are concerned with downwash angles at touchdown, 

E gl * and at nose wheel liftoff, E^, which are proportional to angle of 
attack at these times > and AL^^. The stall velocity V gTaTT at 

liftoff is then found in ft. per sec., and the dynamic pressure at this 
airspeed is Q,—. Subroutine C returns the lift coefficient c T _, 

Ul LIFT LIa) 

and the lift for nose-wheel liftoff is calculated (FL) . 

PL " C LLO C UE S W (V.1.147) 

The moment coefficient at nose-wheel liftoff due to power effects is 
computed in the following lines. This has the form, 

2 D 2 T 
PROP P 

C ■ C — - — E (propeller aircraft) (V. 1.148) 

MPTO TNWLO S w C 0ARW NP 

C ■ f TC/(WC ) (turbojet/turbofan) (V.1.14'0 

for propeller and turbojet aircraft, rsspectively. In these equations, 

the terms C and F are proportional to engine thrust at take-off, 

TNWLO NNWLO 

as were computed earlier in the subroutine. 

Horizontal tail size calculations are initiated following the comment 
card "CALCULATION OF THE HORIZONTAL TAIL SIZE". The tail volume 
coefficients and tail areas are determined for the three horizontal 
sizing criterion. The volume coefficients are functions of previously 
determined stability derivatives and geometric parameters. V and S 

D 1 
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are the volume coefficient and tail area to aatisfy static stability 

and trial criteria. is the neutral point and XCG A p r and XCG^^ the 

allowable c.g. travel limits. V G and are the volume coefficient and 

tail area to meet nose-wheel liftoff criterion. V_„ and S m „ are for the 

GX TX 

stability and liftoff criteria. The most critical of the three criterion 
is determined next and the required horizontal tail area, span, and 
moment arm are computed as 

S HT " V * 3 W CB *V E V (V. 1.150) 

B HT " SQRF (AR HT S HT ) (V. 1.151) 

EL TH " EL H + <Wr ' - 25) CBAR W (V. 1.152) 

V.l.4.2 Vertical Tail Sizing 

The vertical tail is sized to meet a directional stability require- 
ment and in the case of multi-engine aircraft the "engine out" control 

requirement. Basically, the specifications for single-engine airplanes 

do 

state that the directional stability criterion, n , should be negative 

d 

(stable) for any anticipated speed greater than 1.2 times the stalling 
speed. In addition, the yawing moment control, usually the rudder, must 
be powerful enough to counteract the yawing moments encountered in various 
aspects of flight operations such as cross-wind takeoffs or landings. 

These requirements apply also to multi-engine airplanes, with the 
additional requirement that the rudder must be powerful enough to regain 
and maintain straight flight with one engine inoperative at a minimum 
speed not greater than 1.2 times the stalling speed. The minimum speed at 
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which engine-out control can be maintained is known aa the "minimum 

control speed”, V , and the requirement that it be less than 1.2 times 
me 

the stall speed is the major design factor in determining the vertical 
tail and rudder size. Thus, in preliminary design calculations for the 
vertical tail, we are faced with a directional stability problem for 
lignt sing] e-engine airplanes, and a directional control problem for 
multi-engine airplanes. 

The vertical tail sizing computations begin at the comment "VERTICAL 
TAIL SIZING" in subroutine TAIL. The contribution of the various air- 
craft components, wing, fuselage, power, etj., to the directional 
stability of the aircraft are the first calculations. 

The wing contribution to directional stability is in many cases 
negligible. The small contribution of the swept wing is estimated by 

CNP - - 0.00006 DLM. 1/2 (V. 1.153) 

WNG C4 

The fuselage contribution to directional stability is CNP^g, 
is a function of several parameters, as indicated in Figure V.1.28. The 
tail cone and nose cone approximations ~o the fuselage geometry are used 
to derive the equation for projected fuselage side area, S , and the 

b 

fuselage contribution to the directional stability derivative is 

CNP pus - (.96 0/57.3) CS/S/ (EL/B) <D pl /D p2 I 1/2 (V-1 . 154) 

< d f/ d fi |1/3 

where Q is determined from Figure V.1.29, and D and D are fuselage 
B FL F2 

diameters at the quarter and three-quarter body length points, 
respectively. 


V-l 

78 




The wing- fuselage Interference stability increment, which is 
computed next, depends on the wing location, 


DLl - 0.0002 H 
CNP WING 


(V. 1 . 15 r >) 


where “ 1 £or Hi 9 h wing 

■ 0 for Low wing 

The next directional stability parameter computed is the one 

accounting for power effects. This is identified as CNP p and is set to 

zero for turbojet/t.urbofan configurations (N ■ 7). For propeller 

TYE 

aircraft the first item computed is the longitudinal distance front 
propeller disc to center of mass, which is denoted by EL p . This term 
is equal to EL__ for single engine aircraft, since this measures the 
center of gravity aft of the aircraft nose. For twin engine aircraft, 
it is more complex, and is given by 


EL„ - EL_ - <EL_ - .25 + X^) 


CG 


C4W 


(V. 1. 1 5« > 


- eg mrgn * CBAR w + E V 2 ' 


The sideforce stability derivative caused by windmilling propellers 

is implied by the number of blades, B , and the derivative is named 

L 


DC YP 


DP* 


Similarly, the windmilling propeller's contribution to the 


directional stability derivative is CNP po » a function of propeller geometry 
and location, and of wing span and area; i.e.. 


CNP - n * EL DCYP EN / [ 4 . SB) 

PO PROP P DP P W 


(V. 1.157) 
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The full-power value of this derivative is simply estimated by increasing 
the windmilling value by fifty percent. 


C_ _ ® 1*5 C_ 
NPP NPPO 


(V.1.15B) 


The vertical tali moment arm from the aircraft center of gravity is 
computed as 


EL v - el h 


- DEL, 


TVH 


(V. 1.159) 


where EL is the horizontal tail moment arm from the aircraft center of 
H 

gravity and DEL is the distance between the vertical and horizontal 

TVH 

tails aerodynamic centers. 

Another interference correction for directional stability calculations 
is the one that arises from the siuewash or interference from the wing- 
fuselage combination on the vertical tail. This is identified as DL2 CNP 
and is determined from Figure V.1.30. 

The lift-curve slope of the vortical tail, CLA^^, is determined from 
Figure V.1.31 and the "effective" aspect ratio, AR^,^ . of the vertical tail. 

The lateral distance from centerline to engine thrust line is then 
calculated, for use in determining yawing moment due to engine failure. 

This is denoted Y, and it is given in terms of a percent semi-span Y p 
or nacelle diameter DBAR^ and maximum fuselage width SWF. The yawing 
moment coefficient caused by an asymmetric engine arrangement has either 
of two forms, depending on the propulsion type; i.e., for propeller aircraft. 


'NP 


YD C / 
P.?OP TMC 


(.5 


S W B) 


(V. 1.160) 
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EFFECTIVE ASPECT RATIO VERTICAL TAIL 


ar vte 


FIGURE V.1.31 - SLOPE OF LIFT CURVE, VERTICAL TAIL. FROM 
NACA TN 755 "ANALYSIS OF WIND TUNNEL DATA ON 
DIRECTIONAL STABILITY AND CONTROL" H. R. PASS 
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while for turbojet or turbofan aircraft, 


C NP ’ * F NHC ( ' 5 ™°T0 S W B V L> 

The next group of calculations determine the vertical tail areas 
required for the two vertical tail sizing criterion. S QS is the tail 
area required to satisfy the directional stability, CNP ftC , required for 
the aircraft. is the tail area required to meet the engine-out 

requirement for multi-engine aircraft. 


_ ( CNP + CNP_ ttl! + CNP + DL1 „ ♦ DL2 - CNP) S,,B 

S DS - WNG FUS P CNP CNP AC W 


cla lfv ei V et V 


(V. 1.102) 


VT2 


C S B 
NP W 


cla lfv l Vx tau v ei V eta mc 


(V. 1.103) 


where TAU V ■ the vertical tail effectiveness 
DR MAX " maximum rudder deflection, deg 

and all other variables have been previously defined. The most 
critical of the two criterion is determined if it is a multi-engine 
aircrai't and then the span and moment arms are computed for the vertical 
tail. 

(A V S VT ) (V. 1.164) 

+ (X CGAFT " - 25) CBAR W (V.1.10')) 


B - SQRT 
VT 


el tv - EL v 
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V. 2 WEIGHT MODEL USER'S MANUAL 


After the aircraft geanetric characteristics have been established, 
the component weights and balance characteristics are determined by 
straightforward methods. These methods require the use of many descrip- 
tive input parameters which are geometrical, initial or operational in 
nature, and which serve to describe the aircraft mass and its distri- 
bution in certain critical operational configurations. 

The subroutine WGHT , ENGWGT , DLOAD and TAIL, subroutines have 
input-output quantities as defined on the followi g pages. It will be 
noted that WGHT and TAIL have many more inpu* U cput parameters than 
most other subroutines, because of the relatively high level of detail 
required in the weight model analysis. 

Input-output quantities for subroutine WAIT are in Volume IV — 
Propulsion. 
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FIGURE V.2.1 SUBROUTINE DLOAD, INPUT 



CATD 0 to 3, aircraft category indicator 

CBARW wing mean aerodynamic chord, ft. 

CLALPH lift curve slope, per rad 

VMLSFL maximum sea level flight velocity, mi. per hi. 

wos wing loading, lb. per sq. ft. 
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FIGURE V. 2.2 SUBROUTINE DLOAD, OUTPUT 


VARIABLE 

DESCRIPTION 

EMLF 

maneuver load factor, g's 

EMMO 

maximum allowable operating Mach number 

GLF 

gust load factor, g's 

ULF 

ultimate load factor, g's 

VDMIN 

minimum design dive speed, ktu. equivalent airspeed 

VMO 

maximum operating velocity, kts. equivalent airspeed 


V~2 

3 






FIGURE V.2.3 SUBROUTINE ENGWGT, INPUT 



-2 






FIGURE V.2.4 SUBROUTINE ENGWGT, OUTPUT 


DESCRIPTION 


MTGB 


engine specific weight, lb. /lb. thrust or lb./SHP 

engine weight, lb. 

nacelle weight, lb. 

propeller weight, lb. 

pylon weight, lb. 

weight of gear box, lb. 
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FIGURE V.2.5 SUBROUTINE TAIL, INPUT 


VARIABLE 

DESCRIPTION 

ALPHLO 

zero lift angle of attack, deg 

AR 

aspect ratio of wing 

ARHT 

aspect ratio of horizontal tail 

ARVT 

aspect ratio of vertical tail 

B 

wing span, ft. 

BL 

number of propeller blades 

BMLOO 

fineness ratio of tail boom 

CBARW 

mean aerodynamic chord of wing, ft. 

CLMXLD 

maximum lift coefficient, landing configuration 

CLMXTO 

maximum lift coefficient, take-off configuration 

CMFLPL 

pitching moment due to flaps, landing configuration 

CMFLPT 

pitching moment due to flaps, take-off configuration 

CMP LX) 

pitching moment due to power, landing configuration 

CRCLW 

wing chord at fuselage enter line, ft. 

CROOTW 

root chord of wing, ft. 

CXA 

horizontal displacement of main landing gear wheel, 


percent mac 

DBARN 

mean nacelle diameter, ft. 

DELXCG 

center of gravity range, percent mac 

DEMAX 

maximum trailing-edge-up elevator deflection, deg 

DLMC4 

sweep back of wing quarter chord, deg 

DPROP 

propeller diameter, ft. 

DRMAX 

maximum rudder deflection, deg 
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FIGURE V.2.5 SUBROUTINE TAIL INPUT (Continued) 


VARIABLE 

DESCRIPTION 

DWPQCH 

sweep of quarter chord of horizontal tail, deg 

EGMRGN 

engine c.g. relative to leading edge mac, ft. 

ELBM 

length of tail boom, ft. 

ELCG 

distance from aircraft nose to aircraft c.g., ft. 

ELC4W 

position on fuselage centerline of wing quarter chord, ft. 

ELF 

length of fuselage, ft. 

ELODN 

length to diameter ratio of fuselage nose cone 

ELODT 

length to diameter ratio of fuselage tail cone 

ET.PC 

length of pilot compartment, ft. 

ELTH 

distance from wing ac to horizontal tail ac, ft. 

ELTV 

| distance from wing ac to vertical tail ac, ft. 

EMCRU 

cruise Mach number 

ENP 

number of engines 

EYET 

incidence angle of tail, deg 

EYEW 

incidence angle of wing, deg 

GCGMAC 

position of aircraft c.g. in percent mac 

HC 

main cabin height, ft. 

HCK 

cockpit windshield height, ft. 

HN 

• 

nose cone diameter, ft. 

HNCRU 

cruise altitude, ft. 

HWING 

wing location parameter 

KCONFG 

type r>f, fuselage indicator 

KNAC 

nacelle drag accounting indicator 


6 


"i 
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FIGURE V.2.5 SUBROUTINE TAIL INPUT (Continued) 


VARIABLE 

DESCRIPTION 

KWRITE 

output indicator 

LC 

length of fuselage cabin, ft. 

NTYE 

engine type indicator 

RELP 

fuselage-mounted engine location, fraction of fuselage 


length 

RH 

ratio of elevator chord to hori rental tail chord 

RV 

ratio of rudder chord to vertical tail chord 

RVMCS 

ratio of minimum control speed to takeoff stall speed 

SAH 

horizontal tail location on vortical tail 

SLM 

taper ratio of wing 

STATIC 

static margin, fraction mac 

SW 

wing area, sq. ft. 

SWF 

maximum fuselage width, ft. 

TP 

vertical distance c.g. to thrust line, positive downward, 


ft. 

UM 

coefficient of rolling friction of aircraft on dry runway 

WG 

aircraft gross weight, lb. 

WGS 

wing loading, lb. per sq. ft. 

XLN 

1 

length of nacelle, ft. 

XLgDE 

ratio of nacvlle length to diameter 

XMAC 

location of wing leading edge of mac, aft from wing 


centerline leading edge, ft. 
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FIGURE V.2.S SUBROUTINE TAIL INPUT (Concluded) 


VARIABLE DESCRIPTION 

YM^C spanwise location from winy centerline chord to winy 

mac chord 

YP spanwise engine locations on wing/wingspan 
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FIGURE V. 2.6 SUBROUTINE TML, OUTPUT 
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FIGURE V.2.7 SUBROUTINE 


INPUT 


WGHT, 


VARIABLE 

DESCRIPTION 

AR 

wing aspect ratio 

ARHT 

horizontal tail aspect ratio 

AKVT 

vertical tail aspect ratio 

ATMXQC 

maximum tip tank length, divided by wing tip chord 

B 

wing span, ft. 

BHT 

horizontal tail span, ft. 

BVT 

vertical tail span, ft. 

CATD 

1 to 3, specifies structural design criteria 

CBARW 

mean aerodynamic chord of wing, ft. 

CKl-21 

wing scale factors, nominally unity 

CPMRGN 

distance from wing c.g. to aircraft c.g., divided by cbai .. 

DELP 

fuselage cabin pressure differential, lb. per sq. in. 

DELWFC 

incremental weight input to flight control system 

DELWST 

incremental structural weight input 

DLMC4 

wing quarter chord sweep, deg 

r>WPQCH 

horizontal tail quarter chord sweep, deg 

DWPQCV 

vertical tail quarter chord sweep, deg 

EGMRGN 

location nacelle c.g. aft of mac leading edge, divided 
by CBARW 

ELBM 

tail boom length, ft. 

ELF 

fuselage length, ft. 

ELFFC 

cabin length for tail boom fuselage, ft. 

ELINC 

distance from vertical tail LE to horizontal tail LE, ft. 

ELODN 

length to diameter ratio of nose cone 
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FIGURE V.2.7 SUBROUTINE WGHT INPUT (Continued) 



VARIABLE 

DESCRIPTION 

ELPC 

length of pilot's compartment, ft. 

ELRM 

length of pylon attachment for fuselage mounted engines 

EMCRU 

cruise Mach number 

EMLF 

maximum manuver load factor, g's 

EMMO 

maximum operating Mach number 

ENP 

number of engines 

GLF 

maximum gust load factor, g's 

HN 

altitude, ft. 

HCRU 

cruise altitude, ft. 

HPMSLS 

maximum sea level static horsepower 

KNAC 

- 0, nacelle drag hot included in aircraft drag 


■ 1 or 2, nacelle drag included in aircraft drag 

L DC KMX 

maximum allowable tip tank fineness ratio 

NTYE 

1 to 14; type of engine, see ENGWGT 

NTYP 

1 to 16; type of propeller, see ENGWGT 

PAX 

number of passengers 

RELP 

location of propulsion on fuselage, fraction of 


fuselage length 

RELR 

location of c.g. of fuselage and contents, fraction 


of fuselage length 

SAB 

seats abreast in fuselage 

SAH 

horizontal tail location on vertical tail 

SF 

surface area of fuselage 
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FIGURE V.2.7 SUBROUTINE WGHT INPUT (Continued) 


VARIABLE 

DESCRIPTION 

SHT 

horizontal tail area, sq. ft. 

SKB-SKZ 

weight coefficients of aircraft components 

SLM 

taper ratio of wing 

SLMH 

taper ratio of horizontal tail 

SLMV 

taper ratio of vertical tail 

STMRGN 

aircraft c.g. location relative to mac quarter chord 

STRUT 

wing strut attachment point, fraction of semispan 

SVT 

vertical tail area, sq. ft. 

SW 

wing area, sq. ft. , 

SWF 

fuselage width, ft. 

TCHT 

thickness to chord ratio, horizontal tail 

TCR 

thickness to chord ratio, wing root 

TCT 

thickness to chord ratio, wing tip 

TCVT 

thickness to chord ratio, vertical tail 

ULF 

ultimate load factor, g*s 

UWPAX 

weight per passenger, lb. 

VBARH 

volume coefficient, horizontal tail 

VBARV 

volume coefficient, vertical tail 

VDMIN 

• » 

minimum design dive speed, knots EAS 

VMLFSL 

maximum level flight velocity of sea level, mph 

VMO 

maximum allowable operating airspeed, knots EAS 

WENG 

engine weight, lb. 

WFA 

weight of fuel available 
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FIGURE V.2.7 SUBROUTINE VfGHT INPUT (Concluded) 


VARIABLE 


DESCRIPTION 




FIGURE V. 2.8 SUBROUTINE WGHT, OUTPUT 


VARIABLE 1 

DESCRIPTION 


AXIS 

length of tip tank, ft. 

BUT 

horizontal tail span, ft. 

BVT 

• 

vertical tail span, ft. 

BXIS 

diameter of tip tank, ft. 

CBARHT 

mac of horizontal tail, ft. 

CBARVT 

mac of vertical tail, ft. 

CLALPH 

aircraft lift curve slope, per rad 

ELTH , ELTV 

moment arms of horizontal and vortical 

EM 

Mach number 

H 

altitude, ft. 

OWE 

operating weight empty, lb. 

SHT 

horizontal tail area, sq. ft. 

STIP 

wetted tip tank area, sq. ft. 

SVT 

vertical tail area, sq. ft. 

IKPAX 

weight of one passenger, lb 

VFTP 

tip tank fuel volume, cu. ft. 

VFTPMN 

minimum tip tank fuel volume, cu. tfc. 

VFTPMX 

maximum tip tank fuel volume, cu. ft. 

W 

% 

weight of aircraft, lb. 

WB 

fuselage weight, lb. 

WCC 

cockpit controls weight, lb. 

WCFW 

fixed wing controls weight, lb. 

WEP 

total engine weight, lb. 
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FIGURE V = 2.9 INPUT-SUBROUTINE HOPWSZ 


VARIABLE 

AHPSL (I) 

A1IPSLX (I) 

ASWN (I) 
ASWSC (I) 

BMEP 

HPMSLS 

HPQAB 

KSPCHG 

NCYL 

NTYE 

RWH 

SKWDTH 

SKWGT 

XNMAX 



DESCRIPTION 

array of sea level static horsepowers, nonsuperchaigcd 
engines 

array of seal level static horsepowers, supercharged 
engines 

array, normalized specific engine weight, lb. per HP 
array of specific weight increments of engine due to 
supercharger 

brake mean effective pressure, lb. per sq. in. 

maximum horsepower, sea level standard conditions, HP 

horsepower per unit base area, HP per sq. in. 

Z, engine does not have supercharger 
1, engine has supercharger 

number of cylinders in engine (4 to 8) 

engine type (11, 12 or 13) 

ratio of width to height of piston engine 

multiplier for engine width 

multiplier for engine weight 

maximum engine speed, rpcn 



FIGURE V.2.10 OUTPUT-SUBROUTINE HOPWSZ 







FIGURE V.2.11 INPUT-SUBROUTINE RCWS? 


VARIABLE I DESCRIPTION 


Alll’SL (I) 

array o£ sea level static horsepowers, nonsupercharged 


engines 

AH PS LX (I) 

array of sea level static horsepowers, supercharged 


engines 

AROTN (I) 

array of number of rotors 

ASWSC (I) 

array of specific weight increments of engine due to 


supercharger 

ASW75 (I) 

array of specific weights, 1975 technology 

ASW80 (*.) 

array of specific weights, 1980 technology 

GR 

gear ratio 

HPMSLS 

maximum sea level static horsepower 

IDATE 

estimated technology date 

KSPCHG 

0, engine has no supercharger 


1, engine has supercharger 

NTYE 

engine type (14) 

ROTN 

number of rotors 

SKDIAM 

multiplier for engine diameter 

SKWGT 

multiplier for engine weight 

SWSLS 

specific weight of engine (sea- level static), lb. per HP 
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FIGURE V.2.12 OUTPUT-SUBROUTINE RCWSZ 


VARIABLE 


DESCRIPTION 


nacelle surface area, sq. ft. 
gear rati i 

specific weight of engine (sea-level static) , lb. per HP 
length of engine nacelle, ft. 


maximum engine rpm 
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No Input or Output Figures for HOPWSZ or RCWSZ in this volume or the 
propulsion volume. 
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V.3 WEIGHT MODEL PROGRAMMER'S MANUAL 


Aa shown in the flow charts on the following pages , subroutine WGHT 
involves several minor iterative loops and computational sequences 
associated with variati.M. > ' the aircraft geometry. The fuel weight 
(wing and tip tanks) and Uid engine sizing loops are interior to the 
tail area loop, and all of these three major loops must converge before 
the component weights and wing location are fixed. 

The flow charts are given for the aircraft structure weights 
routines WGHT, DLOAD, ENGWGT and TAIL. It may be noted that the input 
to WGHT includes over 70 lengths, angles, and weights which are descrip- 
tive of the aircraft geometry and size. In addition, aerodynamic data 
is required by TAIL, which requires some 85 input parameters, while 
the other two programs are relatively short. Propulsion system weights 
are computed by the routines ENGWGT, HOPWSZ, RCWSZ and WAIT. 
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V.3.1 Structural Weight 


V.3.1.1 Subroutine WGHT, Structural Weight and Wing Poaltioning 
Routine 

This routine controls the calculation of aircraft structural weight 
and the wing positioning calculation. Routines called are DLOAD for 
design load ratios, TAIL for tail sizing based on longitudinal and 
directional stability requirements, CTAER for cruise aerodynamics 
(described in Volume III) , and finally a call tc engwgt passes program 
control to propulsion system sizing. 
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DELWST 

WTIP-0 

WFTP-0 

WFW»0. 

FVELD- 

■0. 

• 

• 

6.687 




■4 or ■11 


NLWGHT=1 
S WPQCH= DWPQCH* . 017453 
SWPQCV=DWPACV* .017453 
RLMC4»DLMC4* . 017453 


CALL DLOAD(VMLi ; SL,WGS,CATD,CLALPH,CBARW,VDMIN,VLF,CLi ; ,i.MI.I ,\MO,I : MMO) 


WFE=61 . 75*SEATS**2. -352. 5 * S I • A I St ; 


WPLL)ES=WPL 

WPLMAX=UWPAX*PAX 


WPLMAX 


=WPLDES 


WPLMAX*WPLDES 


WFWMIN=WFW 
WEP=WENG*ENP 
WPR0P=ENP*WPR0P1 
WPES= (WNAC+WPYLON) *ENP 
WPEI=SKPEI*WEP 

WPSTAR=CK5*WEP+CK7*WPEI+WPROP+WTGB*ENP 


WPES=SKPES*WEP 


WLG=SKLG*WG 

WMG=SKMG*WLG 

CGAMMA= (SW/ 1000. ) **1 . 25* (1 . +VLF*WG/ (200. *SW) ) * (. 01*AR/TCR* , •> ) 

* (1 . +ABS(SIN (RLMC4) )) **. 54* (. 54* (1 . ♦. 2*SLM)* (TCR/ I (' I i **. I 

WC=0. 

WCMIN=0. 

YC=0. 

TSWVLE® (1 . -SLMV) / (1 . ♦SLMV)/ARVT*TAN (SWPQCV) 

CRCLVT»2. *SVT/BVT/ (1 . ♦SLMV) 

XVMAC=BVT*TSWVLE* (1 . * 2 . *SLMV)/3./ (1 . ♦SLMV) 

E LW I NG= ELF - E LTV - CRC LVT* XVMAX+ . 25*CBARVT+CPMRGN*CBARW 
i r 


WMIBU 


FIGURE V.3-1 DETAILED FLOW CHART SUBROUTINE WGHT 











WGHT 


2 


FH*WG*SKY*ELF*BHT*SKTL* (1 . *2. *SL^/(10000007 # (1 . ♦SLMH) ) 

FV-WG*SKZ* (F.LF*B) *BVT* (1 . * 2 . *SLMV)/ (1000000. * (1 . ♦SLMV) ) 

FV*FV/2. 

WHT=350. * (SUT*FH*AL0G10(VDMIN)* (1 . ♦SLMH)/ (200. *ELTH*TCHT*< :n\KIIT) ) '*. 54 


i ai;(iiu ■ ■ axKW ■ u a ; fm u.11 ii : i vkki m t*ai 


/ (200. *ELTV*TCVT*CBARVT) ) ** . 54 


WR=2. * (YP* (WEP*WEPI+WPES*WPROP)+WMG+WCMIN*YC) 


TSWC2=TAN (R1A1C4) - ( 1 . -SLM)/ ( 1 . ♦SLM) /AR 
SWC 2 = AT AN (TSWC 2 ) 

CSWC2=COS(SWC2) 

BSTK*B/CSNC2 

SKNO= 1 . * 2 . 5/SQRT (BSTR) 

SKSTR= 1 . -STRUT ‘STRUT 
SKGEAR= 1 . 0 


SKEPOS= 1 . 0 


KGEAR=0. 95 


‘2 or 3 


SKEPOS= 1 . 05 


WW1 = . 08*WG+WUEI)EV 


SKI:P0S».9:> 


SKEPOS .'in 


E00= (SKSTRMJLE* (WG- . 8*WW1) ) ** . 757 
WW0= (SKWW*SKN0*SKEP0S*SKGEAR*F00*B**1 . 049* (1 . +SLM) 
**. 4)/ ( 1000000. *TCR** . 4*CSWC2**1 . 535) 
WW=WW0+WHLDEV 














WGHT 


3 


i 


k ABS (WW1 -WW) / WWJ 

~u~ 

I wwi=ww I - 




>.005 


■0 


WM=YP/ (YP+. 001)* (WEP*WPEI*WPES*WPROP)+WMG*YMG/ 

(YMG+ . 001 ) ♦WCMIN*YC? (YO . 001) 

WX=WG-WW-WFWMIN-WM 
WBOOM=SKBM*F:LBM*SQRT(XARBM) *WG/1000. 

WB=SKB* ( (WX/ 10000. ) ** . 7* (SF/ 1000. ) *SWF* (ELFFOELRW) ** . 5*AI ltd 1 0 
(VDMIN) * (DELPHI . ) **2 . *ULF** . 3) ** . S08*WBOOM 
WST=CK8*WW+CK9*Wirr^(:K10*WVT+CKll*WB+CK12*WLG*CK14*WPt-S*|)l:l.lvSi*wriP 
QDIV- ( 1 . 15*VI»MIN) **2/391 . 

WSCX=SKFW*SW* *.317* (WG/ 1000. ) ** . 602*ULF** . 525*QDIV** . 345 

WCC=SKCC* ( WG/ 1 000 . ) * * . 4 1 

WCFW=WSCX-WCC 

WSAS=SKSAS 

WFOCKlS*WCOCK18*WCFW*CK19*WSAS*DELWFC 

WFADES= (WG-WPSTAR-WST-WFC-WFE-WFUL-WPLDES)/ (1. ♦CK21*SKFS*(>. <>.S7/FUED) 
/ ( 1 . +CK21 *SKFS*6 .687/FUELD) 

WFSS=CK21* (6. 687/FUELD) *SKFS*WFAD0S 
WP=WPSTAR+WFSS 

WFA=WG-WPSTAR-WST -WFC-WFE-WFUL-WPL-WFSS 

WFAMIN=WG-WPSTAR-WST-WFC-WFE-WFUu-WPLMAX-WFSS 

WFADES=WFA 

WFAMAX=WFADES 


c 


I 


iwr 


J± 


WTEST=WFW 

TC= ( (TCR-SWF/B* (TCR-TCT) ) * (1 . -SWF/B* (1 . -SLM) ) 

♦SLM*TCT) / ( 1 . ♦SLM-SWF/B ( 1 . -SLM) ) 

FVOLW= (SKWF* . 888889*TC* (SW**1 . 5) * (2. *SLM*1 . ) )/ ( (AR** . 5) * ( ( SI.M+ 1 . ) **2. 


)) 


£ 


WFWMX=WFW 

WFXTRA=0. 


c 


WFA 


£ 


> 


< WFW 


Hwfw=wfa ) 


C WFAM1N y — ► I WFWMIN=WFAMIn| 


WFAMIN 


=WFWMIN 


-| WFAMIN=WFWMINj 


WFWM1N=WFW 
WFXTRA= WFAMIN- WFW 
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I 


WGHT 


EXIS-SQRT(AXIS**2.-BXIS**2.)/AXIS 
STIP«2. *(1 . 57*BXIS**2.^1 . 57*AXIS*BXIS 

ARSIN(EXIS)/EXIS)-2. *3. 14*AXMIN*BXMIN/4. 
WTIP-SKWTP*STIP 
WOWTIP^WFTP 
WCMIN-WC 


YC> 1 . 0 


WCMIN«WTIP*WFTRA 


W1=WIT 

A1=WFTP 

IWT*1 


W2=WTTT 

A2=WFTP 

WFTP* (A2*W1 -A1 *W2) / (W1 - A 1 - W2 ♦ \ 2 ) 
IWT*2 


KTI P=0 


LCWING 


XWLG=0. 


NLWGHT=2 


NLWGHT 




WREMN = WG - WHT - WVT - WW - WPSTAR - WF W - WFTP - WT I P 
-WPES- WB00M-CK2 1 *WFSS-XWLG 
ELPROP=RELP*ELF 
ELREMN=ELF*RELR 
CRCLW=2.*SW/B/(1 .+SI.M} 

CRC I.IIT= 2 . * SI IT/ BHT/(1.+SLMH) 

CRC I.VT= 2 . *SVT/ BVT/ ( 1 . ♦SLMV) 

















WGHT 


7 


TSWPLE* ( 1 . -SLM) / ( 1 . ♦SLM) / AR^TAN (RLMC4) 
TSWHLE- ( 1 . -SLMH) / ( 1 . ♦SLMH) /ARHT*TAN (SWPQCH) 
TSWVLE" ( 1 . -SLMV) / ( 1 . ♦SIM) /ARVT^TAN (SWPQCV) 
XMAC«B*TSWPLE* ( 1 . * 2 . *SLM) /6. / ( 1 . ♦SLM) 
XHMAC»BHT*TSWHLE* ( 1 . * 2 . *SLMH)/6. / (1 . ♦SLMH) 
XVMAC«BVT*TSWVLE* ( 1 . * 2 . *SLMV) /6. / ( 1 . ♦SLMV) 


ELC4W=ELF-CRCLVT»XVMAO , 25*CBARVT-ELTV- . 25*CBARW-XMAC» . 2f> *< » I W 

i 


ELPROP 


ELPR0P-ELC4W- . 25*CRCLW*\MAi> 
EGMRGN*CBARW 


CHT=SLM*3. *CBARW/ (2. * (1 . ♦SLM- (SLM/ ( 1 . ♦S^) ) ) ) 
ELTIP*ELC4W- . 25*CRCLW^XMAO . 25*CBARW 
ELT1PE«ELT1P*. 5*AXIS 
ELTPMX=ELC4W+B*TAN(RLMC4)/2. ♦. 75*CHT 


ELTI P=El,TPMX- . 5*AXIS 


KLTIP-O. 


ELCGV=ELC4W- . 25*CRCLW^XMAO . 2S*CBARW*ELTV+ . 0S*CBARVT 

I 

E LCGH= E LCGV - XVMAC - . 25*CBARVT> . 25*CRCLVT>SAH*BVT*TAN (SWPQ(A , 

i 

- . 2 5 *C RC LHT ♦ X1IMAC ♦ . 2S*CBARimELINC 

i* 

YMAC=B* ( 1 . ♦ 2 . *SLMH)/6. / (1 . +SLM) 

i 

YHMAC=BHT* (l . * 2 . *SLMH)/6./ (1 . ♦SLMH) 

■ j? 

h 


jj 


FCGLIM=ELWING- (CPMRGN^ . 25-XCGFWD) *CBARW 
ACGLIM=ELWLNG- (CPMRCN+. S2-XCAFT) *CBARW 
WFUliLM=HPMSLS/2. 


WFUELM= . 2 'WFW 


OWE=WP+WFC+WST+WFE+WFUL 

WREMN 1 = WFUL+WFOWB+WFE+WLG- XWLG 

EMOWE=ELPROP* (WPSTAR+WPES) ♦ELCGH*WHT+ELCGV*WVT+ELWING* (WW* \ iv I .( ;>CK2) 
♦WFSS) ♦WBOOM* (ELF- . 5*ELBM)+ELTIP*WTIP+ELREMN*WKEMN I 
CGOWE=EMOWE/OWE 
NPAX-WPL/200. 

PAS=NPAX 

NROW= (PAX- 1 . )/SAB 
ELCPLT=EL0DN*IlN^ELPC/2 . 

EMPAXT=ELCPLT*170. 

WPAX=PAX*170. 
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<0 


r WFTP 


■0 


(tVFTPD(ILIMIT)^ 

| WFTPU(1 LIMIT) =W1-TP I)( I LIMIT) ♦XFULLD(IL IMIT)) \ 

(wt-mumMin ^ — > ^ 1TP — iwitpihi limit) =wftT*1 


I XIUi:ia)(ILIMir)=WFUi:i>UWFWI)(iLlMrr)-WITPIHlLIMIT)T 


I 


(xRUiUHl LIMIT)) ^ 


WFWI)(I LIMIT) 


> 


<WFWMX 


WFWI)( I LIMIT) =WFW 
Wrm)(ILlMlT)=WITP 

XI' TIE LD ((LIMIT) = WG= WG-OWE - WPAXD ( I L 1 M I T) - WBAG ( I L 1 M I T ) 
-WFWI)( 1 LIMIT) -WFTPD(ILIMIT) 


> 


■0 


-fXFUELIH 1 LIMIT) -0. | 


J(IULL1)( 1 LIMIT) 

WPLIMVPAXIH l LIMIT) ♦WBAG ( 1 LIMIT) 

WFUEL=WFWl)( 1 LIMIT) ♦WFTPD( I LIMIT) ♦ XFUELl)( I LIMIT) 
WT ( I LIMIT) =0WF>WFUI L+WPLl) 


1 


( WT(ILIMIT) - 


-0 


WBAG ( 1 LIMIT) = WBAG ( I LIMIT) - WT ( I LIMIT) ♦WG 
WT ( l LIMIT) =WG 


EML( IL1MIT)=EM0WE*ELWING* (WFW1) (1 LIMIT) ♦XFUELD( I LIMIT) ) ♦LMP \\P( I LIMIT)] 
• +WBAG(I LIMIT) *F.LBAGS+ELTIP*WFTPD(ILIMIT) 

FI LCGI) ( I LIMIT) = EML ( I L I M IT) / WT ( 1 L I M IT) 


88 CONTI Mill: 1 



wiUT!:(o)owi:,a;owi:,wi»AXD(2) ,wpaxd(1) .wpaxd (3) .wbag( 2) , 
lll.BAGS , WBAG ( 1 ) , WBAG (3) ,WFWD(2) ,WFWD(1) ,ELWING,WFWI)(3) . 
WFTPD ( 2 ) , ELTI P .WFTPD ( 1 ) ,WFTPD(3) ,XFUELD(2) ,XFUELl)(l) , 
XFUELt)(3)WT(2) ,ELCGD(2) ,WT(1) .ELCCD(l) ,WT(3) ,ELCGP(3) 
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DELXCG- (ELCGD ( 1 ) -ELCGD (2) ) /CBARW 
ELWING-ELCGD (2) ♦ (CPMRGN+ . 2S-XCGFWD) *CBARW 
ELCP*ELWING-CBARW*CPMRGN 

ELC4W»ELWING-CBARW*CPMRGN-YMAC*TAN(DLMC4*. 017453) 

ELCG=ELWING-CBARW*CPMRGN*STMRGN*CBARW 

GCGAMC= (ELCG-EI-CP+ ■ 25*CBARW) /CBARW 


ELCG= (ELPROP* (WPSTAR+WPES) ♦ (CBARW* (CPMRGN-STMRGN) ) 

* (WW*W! ; W+XWLGh:K 2 1 *WFSS) ♦ELTI P+WFTP) ♦ELCGH* Will 
♦ELCGV*WVT*ELREMN*WREMN+WBOOM* (ELF- . 5*ELBM) ) / 
(WG-WW-WFW-XWLG-CK2 1 *WFSS) 

ELWING=ELCG* (CBARW* (CPMRGN-STMRGN) ) 





WGHT 
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H-HNCRU 

EM=EMCRU 

W=WG 


ICAI.I, CTAl-R 


I GCtlMAC = ( E LCG - E LC P ♦ . 2 S * C B ARW ) / C BAR W 























V.3.X.2 Subroutine DLOAD, Design Load Factors 


This subroutine is celled by WGHT to obtain design speeds and load 
factors in accordance with Federal Aviation Regulations. Design speeds 
and load factors are then employed by WGHT in the structural weight 
analysis as described in Volume V, Part I of this report. 


DLOAD 


r»VMLFSL/l. 15 


CATD 

1 


— 0 

wos 

~T"~ - 

y_>ioo — _ 

j WRITE (6) 

I WOS 


VCAMX= . 9XVMLFKT 
VCCOF=33. 













c 


DLOAD 


3 


€ 


€ 


C. 


i 



RETURN 


c 


4 


C 


e 

preceding page bunk NOT JUMa 

c 




V-3 

20 



V.3.1.3 Subroutine TAIL# Tall Sizing Routing 


Subroutine TAIL determines vertical and horizontal tail aleea on 
the basis of directional and longitudinal stability requirements. These 
requirements are based on Federal Aviation Regulations and current design 
practices as discussed in Volume V, Part I of the present report. 


IVFLG-0 

RHOTO.OC 

DELTVH-EI 

ALWGT-EYI 

DCYPDP-9' 

ARVTE 

CNPAC 


ARVTE* ( 1 . 55+ . 7039*SAH) *ARVT| 


1 99 or XCNPAC 


CNPAC— . 0005*SQRT (WG) /B 


XCNPAC=- . 0005*SQRT (WG) /B 
ZAC=HTG + HWI11G*SWF 
EM=EMCRU 

H»HNCRU 

TPALT (H , GALT , PO , FKALT ,T0 , GO , XKV) ] 

- x — ■ — - ■ a t 

PO=PO*144.0 

CLCR=WGS/(.7*P0*EM*EM) 

♦ 

CLIFT (4 , ALFC ,CLCR, ALPHLO,CLALCR , 999. , AR, DLMC4 ,EM) 

, * mz 

!CLAIjCR+C1.ALCR/57. 3 I 


KENG=5 ! 

A 

ENGINE (THRUST, TSFC, XI ,WF,X2 , PO, TO,EM,H,KENG) 
BLACR=BLANG 

CTCR=THRUST/1 . 4/PO/EM/EM/DPROP/DPROP 
ETACR= ( . 9+ . 1*SAH) » (1 . +8 . »CTCR/3 . 14159) 


ETACR 


ETACR=1.0 


ALP11A=30 . 0 

VSDGE=SQRT (WGS/RHOTO/CLMXTO/. 5) 

EM1=.85*VS0GE/1117. 

HOB=ZAC/B 

* __ 

CLI FT ( 2 , ALPHA , CLMXGE , ALPHLO , CLALLO , HOB , AR , DLMC4 , EMI ) 

— 

CLALLO=CLALLO/57 . 3 

VSIGE=SQRT ( WGS/RHOTO/C tiMXGE/ . b) 

EM=0.85*VSIGE/1117. 

H=0. 0 

KENG=5 

ISEGX=2 


FIGURE V.3.3 DETAILED FLOW CHART, SUBROUTINE TAIL 
































i 


^ KWRITE V‘^~^WR1TE(6,6) VN.QN.ELWnJ 

(j) 

c 


AIL - Pag* 4 



WRITE (6, 5) X4R,X4RN,gftj&,&tAg ^ 


ELH-ELTH- (CBARW*STMRGN ) 

FAR-1./AR-1./(1.+AR**1.7) 

FSLM-(10.-3.*SLM)/7. 

FHT- (1 .-ABS (SAH*BVT) /B) / (2. *ELTH/B) **0.33333 
FLAM-SQRT (COS (DLMC4*. 017453) ) 

DEDAW-4.44* (FAR*FSLM*FHT*FLAM) **1.19 
DEDAW-DEDAW*CLA (AR, DLMC4 , EMCRU) /CLA (AR.DLMC4 ,0.0) 
DCMFUS- (QFUS*SWF**2*ELF) / (SW*CBARW*CLALCP*57 . 3) 
DCMFST- (QFUS*SWF**2*ELF) / (SW*CBARW*CLALLO*57 .3) 
HOB-ZAC/B 

VSTALI.-SQRT (2 . *WGS/CLMXLD/RHOTO) 

EMLD-VSTALL/1117. 


, I — - 

\ CLIFT (6, ALWGL,CLMXLD, ALPHLO, CLALLD, HOB, AR, DLMC4 , EMLD)| 


( CLALLD=CLALLD/57 . 3 

DCMFS L= ( QFUS * SWF * * 2 * ELF ) / ( SW*CBARW*CLALLD« 57 . 3 ) 

, 1 

( \<* y < 2 

^ KELP ^ 


DCMNAC= (QNAC*DBARN**2*ELN) / (SW*CBARW*CLALCR*57 . 3) 
DCMNCT= (QNAC*DBARN**2*ELN)/ (SW*CBARW*CLALLD*57 . 3) 
DCMNCL= (QNAC*DBARN**2*ELN) / (SW*CBARW*CLALLD*57. 3) 



DCMNAC=DCMNAC*ENP 
DCMNCT=DCMNQT*ENP 
DCMNCL=DCMNCL * ENP 

- 

T 


DOMNAC- (VN*QN*ELWN* ( 1 . -DEDAW) ) / ( 28 . 7*ELH*CBARW) 

DCMNAC=DCMNAC*ENP 

DCMNCT=DCMNAC 

DCMNCL=DC!MNAC 


ALPHFT=ALWGT-eyew 

CMFUST=DCMFST*ALPHFT*CLALLO 

CMN ACT= DCMNCT * ALP HFT * C LALLO 

ALPHF-ALWGL-EYEW 

CMFUSL=DCMFS L* ALPHF *CLALLOD 

CMNACL=DCMNCL*ALPHF*CLALLD 

EM=EMCRU 

EMNW=0. 85*VSIGE/1117 . 
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TAIL - Pag* 5 


CLALFT-CLA (ARHT , DWPQCH , EM) /57 . 3 

CLALT0«6. 2832 # COS (DWPQCH*0. 01745) /SQRT (1. - (EM # COS (DWPQCH* .017*15) ) J 

CLALTOCLALTO/57. 3 

CLALTL-CLA (ARHT , DWPQCH , EMLD) /57 . 3 

CLALTT-CLA (ARHT, DOTQCH , EMNW) /57 . 3 


ITRLN (XR, QTAU , RH , TAUH , 7 ) 


ITRLN (XRA,QCA,RH,CHALF, 3) 


ITRLN (XRD, QCD , RH , CHDEL , 4 ) 


CHALH3“CHALF*CLALF'f/CLALTO 
CHDEL3=CHDEL+TAUH* (CHALF3-CHALF) 
EGL-. 5*DEDAW*ALWGL 


WRITE (6, 35) CLALTO,CHALF, CHDEL 


EGT* . 5 * DEDAW * ALWGT 

VSTALL=SQRT ( 2 . *WGS/CLMXTO/RHOTO) 

QUE=.5*RHOTO*(.85*VSTALL)**2 

HOB=ZAC/B 

EM= . 8 5 * VSTALL/ 1117. 


CLIFT ( 2 , ALWGT , CLLO , ALPHLO ,CLALFH , HOB , 


F L=C LLO * QUE * S W 


WRITE (6, 36) QUE,FL 


CMPTO-ENP*CTNWLO*2 . *DPROP**2*TP/SW/CBARW 


DCMCLP=3.*CTCR*DPROP**2*TP/CLCR/SW/(l' 


CMPTO=FNNWLO*TP*CLIO/ (wg*CBARW) 


DCMCLPiO. 















TAIL Page 6 



■j CMFPLT-DFLPTO* (-. 008+0. 003*DFLPTO/40.Tj 


CMFLPL-DFLPLD* (- . 008+0. 003* DFLPLO/40 

i 

CMFPLT-0.0 I 

-p 1 

CMFLPL®0. 0 I 


3 


Tl»XAC-XACF-DCMFUS-DCMNAC-DELXCG-ENP*DCMCLP-S'.ATIC 

T2— (CLALFT/CLALCR) *ETACR* (1 . -DEDAW) * (1 (CHALF3/CHDEL3) *TAUH) 

T3«CLALTL*ETALD*TAUH*DEMAX/CLMXLD 

T4= ( CLALTL *ETALD * (ALWGL-EGL-EYEW+EYET) )/CLMXLD 

T5= ( CMFLPL+CMFUS L+CMNACL+CMPLD ) / (CLMXLD 

V« (T1+T5) / (T2+T3+T4) 

XNO=XAC-DCMFUS-DCMNAC-ENP*DCMCLP-T2*V 

XCGAFT=XNO-STATIC 

XCGFWD=XCGAFT-DELXCG 

ELH=ELTH+ . 25*CBARW-XCGAFT*CBARW 

SST=V*SW*CBARW/ELH 



XCGFS=XCGFWD 

TGl=CMFLPT+CMPTO+CMFUST+CMNACT+CLLO* (XCGFWD-0. 25) 

TG2= ( (Fr.-WG) / (QUE*SW) ) * ( (XG+UM*2CG) /CBARW! 

TG3=CLAETT* (ALWGT-EGT-EYEW+EYET+TAUH*DEMAX) *ETANWL 
VG= (TG1+TG2 ) /TG3 

VGX= (T1+XACF+ (TG2+CMFLPT+CMPTO+CMFUST+CMNACT) /CLLO-0. 25)/ (TG3/CEI/ ' • T2 ) 

STX=VGX*SW*CBARW/ELH 

STO=VG*SW*CBARW/ELH 

VTWST=AMAXl (VGX.VG) 



V=VTEST 

IVFLG=1 

XNO=XAC-DCMFUS-DCMNA'7-ENP*DCMCLP-T2*V 

XCGAFT=XNO-STATIC 

XCGFWD=XCGAFT- DELXCG 

ELH=ELTH+0 . 25*CBARW-XCGAFT*CBARW 

ERR0R=ABS ( (XCGFWD-XCGFS) /XCGFS) 


> 0.01 



ERROR 
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©- 


VERTICAL 



WRITE (6 ,39) Tl ,T2 ,T3 ,T4 ,T5 , V , ELH 
WRITE (6, 40) TGl,TG2,TG3,VG,VGx 


TAIL Pag* 7 



SHT«V*SW*CBARW/ELH 
BHT-SQRT (ARHT*SHT) 
ELTH-ELH*STMRGN*CBARW 
VBARU“SHT # ELTH/SW/CBARW 
CMDEL— CLALTL*V*ETALD*TAUH 
„ - - y - 


WRITE<6,114) 

WRITE (6, 110) 

WRITE (6, 111) ALFC,CLALCR,CLALFT,ETACR,CLCR, DCMFUS , 
DCMNAC , DCMCLP , CTCR 

WRITE (6, 112) ALWGT , CLALLO , CLALTT , ETANWL , EGT , CLLO , 
DCMFST , CMFUST , DCMNCT , CMNACT , CMFLPT , 

CMPTO , CTNWLO 

WRITE (6,113) ALWGL , CLALLD , CLALTL , ETALD , EGL , CLMXLD , 

DCMFSL , CMFUSL, DCMNCL , CMNACL , CMFLPL , CMPED 
WRITE ( 6 , 1 1 5 ) CHALF3 , DEDAW, CHDEL3 , CMDEL , TAUH 


XLEMAC=ELC4W-0. 2C*CRCLW+XMAC 
STANO*X LEMAC+ XNO *CBARW 
STAAFT=XLEMAC+XCGAFT *CBARW 
STAFWD=XLEMAC+XCGFWD*CBARW 


WRITE (6, 100) 

WRITE (6, 101) XNO , STANO , SST 
WRITE (6, 102) STATIC ,STX 
WRITE (6, 103) XCGAFT , STAAFT , STC 
DELXCG , SIIT 
XCGFWD , STAFWD , ELTH 


WRITE (6, 104) 
WRITE (6, 105) 


TAIL 




SIZING 


D* (1 . 0-HWING) *SWF 
CPNWNG— 0 . 00006*DLMC4**0 . 5 
ELNC=ELODN*HN 
ELTC*ELODT*HC 
ELF4+ELF*0 . 25 
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RELP > HELP-^BS (RELP*ELP-ELCGl TAIL Pag* 9 

«<■ ... I 

£ 9797. 


\W 

f BIV (DCYK80,BLACR,BL,BLAAR,BLAR,DCYAR,B, 5,NERR) 1 

I 


DCYKN* ( ( (15.+0. 65*AF) /80.7) -1.0) *0.80+1.0 
DCYPDP*=DCYK80*DCYKN/57. 3 


0. 



WRITE 


(6,54)1 


CNPPO- (3. 141*DPROP**2*ELP i DCYPDP*ENl’) /(4.6 A SW*'bT 
CNPP*1 . 5*CNPPO 


[ ELV=ELH-PELTVH | 


f ITRLN (XDDF , YDEL2 , D/SWF , DL2CNP , 9) ” jj 

I ITRLN (XARVTE,YAV,ARVTE,CLALFV,8) I 

— ■■ I 


c 


Y=YP*B/2 . 0 



WRITE (6,57) DL2CNP , CLALFV, ELV , CNPAC \ 

~ ~ : :_z t 1 



° ■ ■ -■ »< f7=DBAKN/2 . O+SWF/2 . oj — 



Y=0. 0 
CTMC=0 . 0 


NTYE^r-” 1 CNP=(Y«FNMC)/(0.5»R1IOTO»SW«B«VMC»«2) | 


f CNP*(Y»DPROP* f *2*CTMC)/(Q. 5«SW*B)l 


TAUVX=TAUV 


TAUV 


999 


TAUVX . 


ITRLN (XR,QTAU,RV, TAUV ,771 
999., ^ ^ 12.5 



. 0 t ; >WF/2 . 


] 


|tAUV=TAUV* ( 1 ■ - . 0080* '1)KMA~X I .'TsTj 


KWRITE^— WRITE (6, 58) CTMC, Y,CNP,TAUV^ 


CV1=> (CLALFV*ELV*ETAVC)/(SW*B) 

3DS- (CNPWNG+CNPFUS+CNPP+DL1CNP+-DL2CNP-CNPAC) /CV1 

svt=sds 


V- 3 
2UB 












WRITE (6,64) SOS , CNPAC 


TAIL P&ge 10 













I 


1 




t 


I 
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V.3.2 Propulsion System Height 


V.3.2.1 Subroutine ENGWGT , Control Program for Propulsion System 
Weight Calculation 

This routine controls propulsion system weight calculations. Routii.' 
called by ENGWGT are: 

. HOPWSZ - Horizontal opposed piston engine weight and size 
. RCWSZ - Rotary combustion engine weight and size 
. WAIT - Propeller weight (Volume IV) 

. ENGINE - Turbojet and turbofan weight and size (Volume IV) 
Subroutines ENGINE and WAIT are described in Volume IV of this 
report, Propulsion. These descriptions are not repeated here. 


I* 

!? 

:S 

: 5 

!| 

Li 

ij 


9 


f 


t 
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ENGWGT 


WENG=0. 

WNAC=0. 

WPYLON=0. 

WPROP1=0. 

WFAN=0. 

WGB=0. 

WSHR=0. 

NTYE=NTYEX 

NTYP=NTYPX 



SWSLS 


— ^ 1 ,2,3 

SWSLS- 1 . r.* ( 1 . ♦ . 15*KSPCHG) 


| SWSLS= 1 . 5* (1 . ♦ . 15*KSPCHG)| 

f , 



1 SWSLS= . 5 1 




WENG=SWSLS*HPMSLS 
WNAC= VWNAC *ANAC 

WPYI.ON=FPYL* ( (WENC+WNAC) ** . 736) 


1 WPYL0N=0. 


NTYEY=NTYE-10 

KWRITX=KWRITE 



NTYEY 


FIGURE V. 3.4 DETAILED FLOW CHART, SUBROUTINE ENGWGT 






ENGWGT 


2 
































RETURN . 


L 









V.3.2.2 Subroutine H0FWS2, Horizontally Cggowd Piston Engine 


Weight and Sire 

This routine computes weight and size of horizontally opposed piston 
engines. Methods employed are discussed in Section V.l.1.2 of the 
present volume. 



RWL=- . 09*RWII+B 


FIGURE V. 3 . 5 DETAILED FLOW CHART, SUBROUTINE HOPWSZ 
















HOPWSZ 


2 


c 


i 


KSPCHG 


c 


X 


yjo_ 


NCYL 


> 




■4 


c 


x 


c 


NCYL 

T 


> 




NCYL 

f RWL«- . 1S6*RWH+B | 


JbTsH r 

NCYL 

fT 

1 . S3 | 

! 

L- 

=*_ 4 



NCYL 

3 -' ► E 

-Li 1 

J " ■ 

4 




-Lb -i .o J C 

NCYL 

3 ■* -E 

1 • ’i 1 

I 

1 - 

— <+ 



XL=W/RWL 
H-W/RWH 
WN= 1 . 08*W 
HN*1 . 2*H 
XLN* 1 . 2*XL 

AN AC* (WN+2 . 76*HN) * . 8*XLN/ 144 . 
XLNX*XLN 

XLN=XLN/12. 


c 


KRITX 




-►RETURN 


WRITE ( 6 )NTYE, 
KSPCHG, SWSLS, 
HPMSLS , XNMAX , 
BMEP, ANAC.W, 
WN ,II,HN, XL, 
XLNX 


RETURN 


ORIGINAL PAGK 1*3 
< 'F POOR QUALITY 
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V.3.2.3 Subroutine RCWS2, Rotary Combustion Ending Wight and si/ ,< 
This routine computes weight and size of rotary combustion engines. 
Methods employed are discussed in Section V.l.1.3 of the present report. 





DSPRO'l- . 033928* (HPMSLS/ROTN) ** 1 . S 
ENGRPM=91670. * ( (HPMSLS/ROTN) **(- . 5) ) 


CALL B I V (SWSLS , HPSLX , ROTN , AHPSL , AROTN , ASW75 , 8 , 4 , NER ) 


l)SPROT= . 021 1 52* (HPMSLS/ROTN) **1.5 
ENGRPM= 153300. * ( (HPMSLS/ROTN) **(-..' 




CALL BI V (SWSLS, HPSLX, ROTN, AHPSL, AROTN,ASW80, 8, 4, NER) 


RPMPRP=XNMAX*GR 

GR=RPMPRP/ENGRPM 

XNMAX=ENGRPM 


DSPROT 


'^GGINAL Pa GF t© 

>f>' pn^n ® 


P °0R QC t 4lfiy 


|SL=. 173 


|BL=14 . 35 



I m - 1 A 7 A 1 


| DL- 1 O . / O I 


[ BL= 19.18 


BL=21 . 55 



FIGURE V.3.& DETAILED FLOW CHART, SUBROUTINE RCWSZ 




























